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ABSTRACT 


A  ring-reinforced,  orthotropic  circular  cylindrical  shell  subjected 
to  external  hydrostatic  pressure  loading  is  investigated.  The  stresses  and 
deflections  throughout  the  sheil  are  determined  by  a  shell  theory  which  con¬ 
siders  the  combined  and  separate  effects  of  large  rotations,  transverse  shear 

1  I 

deformation,  initial  deflections  and  Fiugge  type  thickness  terms.  The  ring 
deformations  are  described  by  both  a  deep  ring  theory  and  a  plane  strain 
analysis.  The  results  are  used  in  the  development  of  pertinent  design  for¬ 
mulae.  Numerical  results  applicable  to  a  typical  glass-reinforced-plast ic 
shell  indicate  that  the  nonlinear  effects  may  be  accurately  provided  for  by 
using  a  perturbation  solution-  Hill's  criterion  for  yielding  of  an  orthotropic, 
material  as  well  as  an  analysis  to  approximate  the  actual  stresses  in  the  con¬ 
stituent  materials  of  a  nonhomogeneoi i  shell  are  also  applied. 
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parameter  measuring  higher  order  effects  in 
governing  differential  equation,  Eq.  (2.4) 

arbitrary  constants  of  integration  in  Lame 
analysis  of  ring,  Eqs-  (3-10) 

constants  in  elastic  law  relating  stress  re¬ 
sultants  to  strains  and  curvature  change,  Eqs.  (2.13) 

cross-sectional  area  of  ring  frame,  Eq.  (3-3) 

parameter  measuring  coupling  of  beam  column 
effect  with  transverse  shear  effect,  Eq-  (2.4) 

arbitrary  constants  of  integration  of  solution 
to  governing  differential  equation  of  deflection, 

Eq.  (4-7) 

axial  flexural  rigidity  of  shell,  Eq -  (2.17) 
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modified  flexural  rigidity  (includes  Flugge 
and  rtonhomogeneous  effects),  Eq.  (A- 2.8) 

moduli  of  elasticity  in  axial  and  circumferential 
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modified  elastic  constants,  Eq.  (A  28) 
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transverse  shear  modulus 

modified  shear  modulus,  Eq.  (A - 28) 

constants  measuring  ratio  of  ring  deflection  to 
interaction  pressure  corresponding  to  ring  theory, 
isotropic  and  orthotropic  Lame  analyses,  respectively 

length  of  unsupported  shell  betwe.en  rings.  Fig.  I 

center  to  center  distance  between  rings,  Fig.  1 

effective  width  of  shell  plating,  Eq-  (6 - 1 0) 

axial  bending  moment  per  unit  circumferential  length 

axial  force  per  unit  circumferential  length 

circumferential  force  per  unit  axial  length 

transverse  shear  force  per  unit  circumferential 
length 

Q.  evaluated  at  the  frame  (i-e-,  at  x  =  ±  L/2) 

equivalent  radial  line  load  (per  unit  circum¬ 
ferential  length)  on  ring,  Eq.  (6.2) 


equivalent  line  load  on  section  of  shell  in 
contact  with  frame,  Eq.  (6-6) 
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radius  of  datum  circle  in  ring,  Fig.  1 
radial  coordinate,  Fig.  1 

radii  of  outer  and  inner  ■  ing  surfaces.,  rerpert  :  vc  1  v,  Fig.l 

radii  of  interior  surfaces  of  outer  and  inner 
flanges  of  ring,  respectively,  Fig.  1 

radius  of  surface  of  contact  between  ring  and 
shell,  Fig.  1 

nond imens iona I  shell  deflection  parameter,  Eqs.  (5-3) 

W  at  midbay  and  frame,  respectively 

axial,  circumferential,  radial  and  transverse 
shear  yield  stresses,  respectively 

radial  ring  coordinate  measured  positive  inward 
from  ring  datum  circle.  Fig.  1 

mean  shell  radius.  Fig.  1 

constants  in  elastic  law  relating  stress  resul¬ 
tants  to  strains  ant  curvature  change 

faying  width  of  ring,  Fig.  2 

width  of  ring,  Eq.  (3-3) 


width  of  ring  web 


widths  of  outer  and  inner  flanges  of  ring, 
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nond imens iona 1  geometric  ring  parameter,  Eq.  (J.4) 

depth  of  the  ring  cross  section,  Fig.  2 

plastic  potential  function  of  the  stresses,  Eq.  (7-1) 

characteristic  functions  of  shell  differential 
equation.  Eqs .  (4.8) 

characteristic  function  evaluated  at  x  =  L/2 
shell  wall  thickness,  Fig.  1 

rad;al  distances  from  dattsn  circle  to  outer  and 
inner  surfaces  of  ring,  respectively,  Fig.  1 


P 

p  ,  p',  p.,  p! 

ro  o  t  ri 

P 


P 


cr 


externa!  hydrostatic  pressure 

pressures  acting  on  ring  surfaces.  Fig.  2 

axisymmetric  buckling  pressure  of  an  unstiffened 
orthotropic  circular  cylindrical  shell  under 
axial  pressure,  Eq.  (6.15) 

axisyrmietric  buckling  load  of  ring  supported 
orthotropic  circular  cylindrical  shell  under 
hylnwstatic  pressure  (6. 16) 


interaction  pressure  between  ring  and  shell, 
Eq-  (6.4) 
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pressures  required  for  shell  to  yield  at  the 
position  (x,z)  predicted  by  Hill's  criterion 
and  maximu.  stress  criterion,  respect  1 . „ 1 , , 

Eq-  (7-2) 

approximate  value  of  p^(x,z),  Eq .  (7-4) 

pressure  at  which  yielding  has  penetrated  through 
half  the  shell  thickness,  Eq.  (7-5) 

lengths  of  infinitesimal  line  element  in  shell 
wall  prior  to  and  after  loading,  respectively 

thicknesses  of  outer  and  inner  flanges  of  ring, 
respectively 

axial  displacement  of  shell  median  surface 

axial  displacement  of  point  in  sheii  wail 

inward  radial  deflection  of  shell  median 
surface  measured  from  median  surface  prior 
to  loading 

initial  deflection  of  median  surface 


inward  radial  displacement  of  ring  outer 
surface^  Fig*  1 

inward  radial  displacement  of  length  of 
shell  in  contact  with  ring 

particular  solution  of  governing  differential 
equation 

ax’al  and  inward  radial  coordinates  measured 
from  the  median  surface  of  the  deformed  shell 
at  midbay;  Fig.  1 

value  of  x  at  position  of  minimum  bending 
stress,  Eq.  (6*9) 

parameter 

inward  initial  deflection  at  midbay 

initial  deflection  for  which  shell  particular 
solution  vanishes,  Eq.  (4.2b) 

initial  deflection  for  which  the  shell  deforms 
without  bending,  Eq  {4.10c) 


initial  deflection  for  which  interaction 
pressure  vanishes,  Eq .  (6.4) 

constants  depending  on  yield  stresses  and 
Poisson  ratios  of  shell  material 


.  nondimensional  shell  curvature  parameter, 

Eq-  (5-5) 
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values  of  A  at  midbay  and  at  frame,  respectively 

ring-shell  interaction  parameters 

parameter  measuring  higher  order  effects  in 
governing  differential  equation 

values  of  i  at  buckling  load 
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parameter  measuring  Flugge  effect,  Eq.  (2,4) 

nondimensionai  ring-shell  interaction  parameters 
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measure  of  the  beam-column  effect,  Eq.  (2.4) 
value  of  y  fof'Ahich  shell  buckles 
nondimensionai  load  parameter,  Eq .  (2.4) 
transverse  shear  strain 

=  h/2a 

internal  or  external  ring  indicator 

axial  and  circumferential  strain  components, 
respect i ve ly 

axial  and  circumferential  median  surface  strain 
components,  respectively,  Eq.  (2.12) 
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contact  with  the  ring 
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radial  ring  strain 

inward  radial  and  axial  coordinates,  measured 
from  midbay  at  the  median  surface  of  a  perfect 
cylindrical  shell,  to  points  on  the  shell  before 
Oseqo  subscript)  and  aft,er  (one  subscript)  de¬ 
formation,  Eqs.  (A.I),  (A .  1 ) 

real  and  imaginary  components  of  roots  of 
characteristic  equation,  Eq.  (4-5) 
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roots  of  characteristic  equation 

nondi mens ionai  ring  material  pa> ameter,  Eq.  (4.4) 

nondimens i onal  ring  shell  interaction  parameter 

Poisson  ratios  of  shell  material 

modified  Poisson  ratio,  Eq.  (2.J7) 

ring  Poisson  ratios 

nondimens ionai  measure  of  the  transverse  shear 
effect 


axial  and  circumferential  shell  normal  stresses 

axial  and  circumferential  bending  stresses 

radial  and  circumferential  normal  stresses  in 
r'itlg,  despect'iveiy 

shell  circumferential  membrane  stress  due  to 
axial  bending 

transverse  shear  stress 


circumferential  cylindrical  coordinate 


axial  rotation  (measured  with  respect  to 
initially  deformed  cylinder) 


I.  INTRODUCTION 


I. 


Attempts  to  design  submersible  vehicles  capable  of  extended  lateral 
excursions  at  great  depths  have  led  to  extensive  research  on  ring-reinforced 
cylindrical  shells  constructed  of  filament  wound  composites.  This  material  is 
relatively  eat/  to  fabricate  and  possesses  a  much  higher  strength-to-wei ght 
ratio  than  the  commonly  used  metals.  The  shells  (and  ring-supports)  are  formed 
by  winding  strong  glass  fibers  about  a  cylindrical  mandrel  and  either  impreg¬ 
nating  or  preimpregnating  them  with  an  epoxy  resin. 

For  many  applications  the  shell  material  as  described  above  may  be 
considered  as  orthotropic  and  non homogeneous  through  the  thickness-  For  deep- 
submergence  vehicles  the  shell  thickness  must  be  relatively  large  in  order  to 
withstand  severe  pressures.  Observed  shell  failures  indicate  that  transverse 
shear  may  be  a  significant  factor  of  the  collapse  mechanism  (Ref.  1).  Further¬ 
more,  built-in  outward  initial  deflections  may  enable  the  designer  to  utilize 
the  "beam-column"  effect  to  strengthen  the  shell.  Thus,  the  analysis  presented 
here  considers  the  separate  and  combined  influence  of  these  effects  on  the 
stress  state  of  the  shell. 

Numerical  results  demonstrate  that  the  more  accurate  calculation 
of  the  cross-sect iona 1  area  of  the  shell  (which  increases  the  axial  stress 
resultant)  is  the  only  significant  thick  shell  correction.  It  is  found  that 
the  stresses  do  change  when  transverse  shear  is  considered.  However,  the  com¬ 
parison  made  wi th  the  three-dimensional  elasticity  analysis  developed  in  Ref.  2 
indicates  that  no  improvement  is  obtained  by  permitting  transverse  shear  de¬ 


formations. 


2. 


Because  of  the  low  radial  yield  stresses  of  the  composite,  yielding 
(as  described  by  Hill's  orthotropic  criterion)  can  begin  well  in  advance  of 
experimentally  obtained  collapse  loads.  The  elastic  region  is  thus  confined 
to  relatively  low  pressures-  This,  in  turn,  keeps  the  deflection  small  enough 
to  avoid  any  significant  beam-column  effects.  Initial  deflections,  on  the 
other  hand  may  be  large  enough  to  cause  substantial  beam-column  effects  even 
at  relatively  low  pressures-  Thus,  initial  deflections  could  affect  the 
stresses  and  hence  change  the  value  of  the  load  for  which  yielding  begins- 
To  obtain  accurate  estimates  of  the  stresses  in  the  constituent 
materials  it  is  necessary  to  consider  the  nonhomogeneous  nature  of  the 
structure.  When  the  stresses  are  thus  obtained  and  utilized  in  a  maximum 
stress  yield  criterion,  the  yield  pressures  obtained  are  reasonably  close 
to  the  experimental  collapse  pressure. 

The  deformation  of  the  reinforcing  frames  are,  in  general,  adequately 
described  by  ring  theory  (which  permits  only  circumferential  stress).  However, 
the  frames  of  interest  here  are  constructed  of  an  orthotropic  material  for  which 
radial  deformations  may  become  significant.  For  this  reason  a  more  accurate 
plane  strain  analysis  is  performed,  in  addition  to  the  ring  analysis,  and  a 
comparison  is  made  for  interna!  frames  of  rectangular  cross  section.  Excel¬ 
lent  agreement  is  found  for  a  particular  ring  considered  in  this  work. 

Designers  of  present  day  submersibles  are  familiar  with  the  for¬ 
mulae  first  developed  for  thin,  isotropic  metal  shells  by  van  Sanden  and 

p  i 

Gunther  (Ref.  3)>  which  were  later  modified  by  Viterbo  tRef.  4),  and  finally 
by  Pulos  and  Salerno  (Ref.  5)-  The  results  presented  here  represent  an  ex- 


3- 

tension  of  the  work  of  Ref-  5  in  th?t  the  formulae  obtained  therein  for  use 
of  metal  shells  are  modified  to  include  more  complicated  filament  wound  com¬ 
posite  she'ls.  These  include  solutions  for  various  stresses,  points  of 
minimum  bending,  ring-shell  interact  ■  or.  pressure,  axisymmetric  buckling 
loads,  effective  width  of  shell  plating,  stress  in  the  reinforcing  rings, 
yield  pressure  and  the  location  at  iihich  yielding  begins. 


4. 


2.  Governing  Shell  Equations 

The  differential  equations  of  equilibrium  of  a  circular  cylindrical, 
orthotropic,  rinq~supporred  shell  subjected  to  hydrostatic  pressure  loading 
are  found  to  be  [see  Fig.  1  and  Appendix  A,-  Eqs.  (A. 20),  (A.3l)  and  (A. 3?)] 

",****  +  Me/L)2(A/B)wiXX  ♦  M9/l)4w/b  -  ('•re/»w0)WXX 

-  [4y(6/L)2/Bjwo^xx  +  8(0/L)27ll-(v^/2)(l+6)]/aB(l+6)  (2.1) 

«  =  [Sw^xx-4|  (©/L) 2w^8r^[  l-(v^/2)(l+6)J/a(l+6)-4yewo  xxj/(l+2a|)  (2.2) 

The  associatedcond i t ions  at  the  shell  boundaries  x  =  ±  L/2  are 

w  =  ws  =  wft  (2-  3a) 

a>  0  (2. 3b) 

In  Eqs.  (2.1)  to  (2.3)  w  is  the  inward  deflection  of  the  shell  median  surface 
from  the  initially  deformed  shape  and  o>  is  the  additional  rotation  of  the 
median  surface  normal  about  a  circumferential  tangent;  w  is  the  inward 
radial  ring  displacement  at  the  contact  surface  between  the  ring  and  the 
shell  and  wg  is  the  inward  deflection  of  the  section  of  shell  in  contact  with 
the  ring. 
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3-  RING-SHELL  INTERACTION 

The  rings  considered  are  assumed  to  be  constructed  of  a  linearly 
elastic  homogeneous  material  and  to  have  cross  sections  of  the  generalized 
"I"  type..  If  a  ring  is  assumed  to  resist  its  loading  with  only  a  circum¬ 
ferential  stress,  the  corresponding  strain  is 

«s  =  “  wR/(*  '  Z)  (2-0 

where  wR  is  the  radial  displacement  (positive  inward);  R  is  the  radius  of 
a  datum  circle.  Z  is  the  radial  coordinate  measured  positive  inward  from 
the  datum  circle.  The  ci rcumferentia  1  stress  is  given  by 

-  £sf's  -  -  -  z> 

E  ^  is  the  modulus  of  elasticity  in  the  circumferential  direction-  The  cor¬ 
responding  hoop  force  is 

h  i 

Fs  =  J  osb(Z)dZ  =  -  E3fAfc'(wR/R)  (3.3) 

-h 

o 

where  b(Z)  is  the  width  of  the  ring  at  Z;  hQ  and  h.  are  the  radial  distances 
from  the  datum  circle  to  the  outer  and  inner  ring  surfaces,  respectively, 
and  A^  is  the  ring  cross-sectional  area.  The  nondimensional  parameter  c'  is 
defined  by 

h. 

=  f-J*  ‘  b(Z)/(R  -  Z)dZ 

f  -h0 


C 


(3.4) 
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and  depends  only  upon  the  geometry  of  the  ring  cross  section.  Appendix  B, 
taken  from  Ref.  .7.,  evaluates  c1  for  different  types  of  "j"  ring  cross 
sect  ions • 

The  hoop  force  acting  on  the  ring  is  easily  obtained  by  cutting 
the  ring  in  half  and  finding  Fg  necessary  to  maintain  equilibrium  (Fig.  2b). 
This  leads  to 
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(3 -5a) 


where 


b' 
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(3 -5b) 


b  J  ■  b .  -  b 
i  1  w 


(3- 5c) 


is  the  web  width,  bQ  and  b.  are  the  outer  and  inner  flange  widths,  res¬ 
pectively;  and  t  and  t.  are  the  thicknesses  of  the  outer  and  inner  flanges, 
respectively.  po,  p^,  p.,  p!  denote  pressures;  the  subscripts  o  and  i  de¬ 
note  the  outer  and  inner  flanges,  respectively,  while  the  prime  superscript 
denotes  the  surface  of  the  flange  nearest  to  the  datum  circle  (see  Fig.  ?) . 

If  the  frame  acts  on  the  inner  shell  surface. 


(3- 6a) 


Fo  = 


(3.6b) 
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b  =  b 
c  o 


(3.6c) 


If  the  frame  acts  on  the  outer  shell  surface, 


p  1  =  p 1  =  p  =  p  • 

r  i  ro  ro 


p.  =  p 
I  c 


b  =  b. 
c  i 


(3.7a) 

(3.7b) 

(3.7c) 


Pc  is  the  interaction  or  contact  pressure  between  the  ring  and 
the  shell.  b£  is  the  faying  width  of  the  contact  surface  between  the  ring 
and  the  shel 1 . 

If  Fs  is  eliminated  between  Eqs.  (3-5a)  and  (3-3),  a  relation  be¬ 
tween  the  radial  ring  displacement  and  pressure  loads  is  found 


WR  =  “  ^cKK  +  piRibi  ‘  PoRobo'pibi’Ri]/(EsfAfc,) 


(3.8a) 


where 


R‘  =  R  -  t 
o  o  o 


(3-8b) 


*1 


(3.8c) 


and  Rq  and  Rj  are  the  radii  to  the  outsr  and  inner  ring  surfaces,  respectively. 


10. 


A  more  accurate  description  of  an.  orthotropic,  ring  of  rectangular 
cross  section  may  be  obtained  by  permitting  radial  strains  in  addition  to 
hoop  strain  via  a  Lame  plane  strain  analysis.  The  field  equations  necessary 
to  describe  such  ring  deformations  are 


CR  "  "  WR,R 

(3*  9a) 

«.  =  -  VR 

0-9b) 

eR  =  (°r  “  VRsfas)/ERf 

0-9c) 

€s  =  1  (°s  *  VsRfaR)/Esf 

0-9d) 

^R°R),R  '  °s  =  0 

0-9e) 

R  is  the  radial  coordinate  measured  from  the  center  line  of  the  ring,  b  ,, 

Rr 

E  c,  v_  ,  v,  are  the  elastic  constants  of  the  ring.  The  solution  is 
Sf  Rsf  SRF 

cR  =  ERfl-  A(X  +  vsR^)RX_'  +  -  vsR  )R“X“,]/(1  -  vf2)  0.10a) 

a£  =  E.fl-  A(x.vft$^  +  1)R-X-'  +  B(XVr<.  -  DR^”1 3/(1  -  v/)  0.10b) 

wR  =  ARX  +  BR_X  0. 10c) 


where,  for  an  internal  frame 


1  1  ■ 

A  =  [pc(|-Vf2)R01-XJ/tERf(^vsRf)t  t-(R|/Re)2^)j 

(3.10d) 

B  -  lpc(l-vf2)S0,'fXJ/U((f(»--vs|,  )[(VKl)2k-'D 

(3- I0e) 

^  ‘  £Sf/ERf 

(3.iof) 

vf2  '  vsfi,vl(s< 

(3  lOg) 

Thus,  for  an  interna  1  frame  the  relation  between  interaction  pressure 
and  the  radia'  displacement  of  the  frame  outer  surface  is  expressed  by 

«„  fVtV'1*2  < 1  -vf!)  [  <X-Vs«f <W'’,R{>  < 1 

:  i (x+vsR  ) (X-V^R  )[ l-(l-d/R0)2xJ]  «  (3.11) 

where  d  is  the  depth  of  the  ring  cross  section. 

For  an  isotropic  materia! 

(EsfwR/pcRo)  =  l(!-Vf)  +  (!+vf)(l-d/Ro)2]/U-(l-d/Ro)2J  =  Kj  (3-12) 

For  an  internal;  rectangular  cross  section  frame,  ring  theory  gives 
'see  Eq  .  (3--  8)  ] 


(Esfw/pcRo)  =  -  (2-d/Ro)/((d/Ro)[2+(f/6)(d/Ro)2(l-d/Ro)2l]  *  K) 


(3-13) 


12. 


A  comparison  of  the  two  solutions  may  be  obtained  from  Eqs.  (3-12)  and  (3- i 3) 

1  -  (K2/K,)  w  v(jf)  +  Of  (~)2  J  (3.Ut) 

O 

Thus,  for  thin  isotropic  rings  the  two  solutions  are  in  very  good 
agreement.  Furthermore,  when  the  orthotropic  Lame  solution,  is  compared 
with  the  ring  solution,  the  agreement  is  even  better  for  the  particular  shell 
examined  in  this  work. 

The  radial  ring  displacement  and  shell  deflection  must  be  matched 
at  the  shell  boundaries.  Because  the  shell  theory  used  here  cannot  describe 
the  radial  changes  in  stresses  that  occur  in  the  shell  region  in  contact  with 
the  ring,  that  length  of  the  shell  is  treated  separately  with  a  simple  analysis- 
Since  the  ring  and  shell  are  bonded  together,  and  classical  ring  theory  does 
not  provide  for  distortion  of  the  ring  flanges,  the  assumption  is  made  that 
the  portion  of  the  shell  in  contact  with  the  ring  does  not  bend.  With  this 
assumption,  equilibrium  of  radial  forces  leads  to  (see  Fig.  2c) 

bcNs  -  2$fa  +  [ ( 1+6) p6 ”  +  (1+5-26  )pc]abc  =  0  (3-15) 

where 

,1  for  internal  frame 

£'  =  { 

0  for  external  frame 


(continued  on  next  page) 
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>c  •  {>“ 


b  for  internal  frame 
o 

for  external  frame 


and  from  £q.  (A. 27b)  and  the  above  assumption 


w 

-  (N  -  v"  N  )  /hE 
a  s  qix  x <p 


Combining  Eqs.  (3-15)  and  (3-17)  to  eliminate  Ns  yields 


ws  =  (aVhEp[  (1+6)P6'C  +  P<_ ( 1 +£ -26")  -  (2<if/bc)  +  v*  (Nx/a)] 


The  matching  conditions  are  that 


wR  =  ws 


ws  =  wf 


Equation  (3 - 1 9a)  is  used  to  satisfy  the  boundary  condition 
(A. 20b).  Eq.  (3-19*),  combined  with  Eqs.  (3-18)  and  (3-13),  leaves  an 
pression  for  the  interaction  pressure  in  terms  of  quantities  which  can 
determined  from  the  shell  solution 


Pc  =  (ac/PcH'(2Vt>c)  +  V(Nx/a)  +  Pl0+6)fi*+A*]j 


(3-16) 

(3-17) 

(3. 18) 

(3-l9a) 

(3.19b) 

Eq. 

ex- 

be 


where 


14. 


X*  =  (Pc/ac)(l-6V')[(Af/bcRc)  -  1] 

\  -  yq 

“c  ■ 


for  internal  frame 
for  external  frame 


P*  =  (26*- 1) (a.  +  3c)  -  0.6 


(3  20) 


The  term  in  Eq.  (3-17)  corresponds  to  the  "Viterbo  effect" 
(Ref.  4)  in  that  it  accounts  for  the  two-dimensional  nature  of  the  stresses 
in  the  straight  section  of  the  shell  adjacent  to  the  ring. 
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4.  SOLUTION  OF  GOVERNING  SHELL  EQUATIONS 

The  solution  of  £qs.  (2.1)  and  (2.2)  subject  to  boundary  condition 
Eqs .  (2.3)  can  now  be  obtained.  Lunchicfc  and  Short;  Ref.  (8)  suggested  that 
the  initial  deflections  encountered  in  isotropic  metal  shells  could  be  approx 
mated  by  a  simple  parabolic  function  which  would  ensure  constant  axial  curva¬ 
ture.  This  form  is  assumed  here;  thus 

wo  =  a[ 1  -  4(x/L)2]  (4.1) 

where  A  is  a  constant  equal  to  the  largest  value  of  wq  (at  midbay) .  With 
the  type  of  initial  deflection  described  in  Eq.  (4.1),  a  particular  solution 
to  Eq .  (2.1)  is 

wp  =  8(7A  -  7Ap)/02  (4.2a) 


where 


Ap  =  -  L2[l  -  (v*/2)(l  +  6)]/4a(l  +  6)  (4.2b) 


If  the  complementary  solution  of  Eq.  (2.J)  is  assumed  to  have  the 

form 


w 


Xx 


(4-3) 


J6. 


The  characteristic  equation  becomes 


\k  +  Me/L)2(A/BU2  +  4(9/L)4/B  =  0 


(k.k) 


with  the  four  roots 


X1,2,3,U  =  ± 


or 


(^■5a) 


X)j2,3,4  =  ±  2^e/LHTl  ±  'V 


(^•5b) 


where 


T) 


1,2 


(«» -5c) 


if  =  A/'/B 


(4.5d) 
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Numerical  results  demonstrate  that  for  the  range  of  £  and  a 

of  interest  in  this  report,  the  shell  reaches  a  yield  condition  while 

2 

7,  £  and  a  are  in  the  region  defined  by  ♦  <1.  Moreover,  in  this  re¬ 

gion  the  following  special  cases  can  be  obtained; 

1.  Negiect  beam-column  effect  by  setting  7=0 
2-  Neglect  transverse  shear  effect  by  setting  £  =  0 

1  I 

Neglect  Flugge  type  corrections  by  setting  a  =  0  and  6  =  0. 

Therefore,  the  only  solution  considered  in  this  work  will  satisfy  the  in- 
2 

equalities  0  c  i|r  <1  and  B  >  0.  The  final  form  of  the  solution  depends  on 
whether  : 


B  -  0 

< 


(4.6a) 


*  =  0 


(4.6b) 


0  <  i|r  <1 


(4.6c) 


where 


Thus,  the  solution  of  Eq.  (2.1);  can  be  written  in  the  real 


w(x)  =  wp[ 1  +  I  c  j  g  j (x) ] 


gj  (x).  *  s  inh  (2T)|0x/L)cos (2T|20x/L) 


(continued  on  next  page) 


g-(x)  ■  *  inh  (2Tl|®x/L)s  in  (2Tl_ex/L) 


g^(x)  *  cosh(2Tij0x/L)sin  (ITl^Sx/L) 

g^(x)  =  cosh  (2Tlj0x/L) cos  ^T^tSx/L)  (4.8) 

Cj  represents  four  arbitrary  constants  of  integration.  Symmetry  demands 
that  (see  Fig.  1) 

c,  =  c^  =  0  (4. 9) 

The  constants  c^  and  are  determined  by  using  the  boundary 
conditions  defined  by  Eqs.  (2-3)  in  conjunction  with  Eqs .  (3- 18),  (3-20), 
(A. 34),  (A-35),  (4.2),  (4-7)  (48),  and  (4-9).  Hence, 

c4  =  -  A "(7b  -  7A j ) / ( 7A  "  7Ap)  (4.10a) 


where 


A*  =  P*B(^/tBP*+G1(1-Pc)[p*+(l+6-2B*)ac]/Pc)  (4.10b) 

A,  -  Ap+L2([ l+( 1+6-26*)  (Oc/p*) ][B*-(v  /2)  (1+6) ] 

+  a  X'7[p“(l+B) ]j/4a(l+B)  (4.10c) 


where 


G,  =  l(Ve)(9J0+g;J0)]/[(B1/T)l)(g50g20+g)0g40)-rvB2/n2)(g50g40-g)0g20)j  (4.12a) 

G2  =  t  (B|/n|)g,0+(B2/Tl2^930''/^  ^30920+910940^  +  ^B2/T12^  ^93094o"9  1 0920^  ^ 

(4. 12b) 


l  ( 1  /n2)  ^4o95o'9 i 0920)  (b  1  ■ ^iO - 1 1 /T1 1 )  (9}0920+9 1 0g40)  }/L  (Bj/n , )  (g3o92o 

+  9io940)  +  '  9i0920)] 

(4. 12c) 

L(B,/Tl2^930  "  ^B2/,T*1^9I0^/^B|'/91^930920+910940^ 

+  (e2/92)(930940"9l0920)] 

(4. 1 2d) 

B,  =  1  -  2 gl  21  -  2a  -  (1-2*)  ] 

(4. 12e) 

B2  =  1  -  2§  [  2 £,  -  2a  +  (1+2*)] 

(4.l2.f) 

910’920’930’940  =  91(L/2),92(L/2),g3(L/2),g4(L/2) 

(4 .  12g) 

20. 


The  junction  G^  is  presented  orematurely  here  in  order  to  list 
the  complete  set  of  Gj,  which  is  used  subsequently.  The  functions  are 
similar  to  ihose  plotted  as  functions  of  y  and  6  for  the  isotropic  case 
[see  Ref.  5  and  App -  D,  Eqs.  (D12)].  where  the  effects  of  transverse  shear 

I  I 

and  Flugge  terms  were  not  considered  (i.e.,  I  =  6  =  a  -  0).  When  these 
effects  are  accounted  for,  these  functional  relationships  change  in  that 
G  j  =  G  j  (7,  0,  l,  L,  Of)  ■ 

It  is  interesting  to  note  that  in  their  work  on  isotropic  shells 
Lunchick  and  Short,  Ref.  (8),  neglected  transverse  shear  effects  and  found 
a  "critical"  initial  deflection  parameter  such  that  the  shell  deforms  with¬ 
out  bending  (i-e.,  cn  =  0) .  This  parameter,  modified  here  to  include 
,x 

orthotropic  effects  and  thick-shell  terms,  which  results  from  consistently 

I  I 

maintaining  Flugge  accuracy,  is  now  given  by  A(  [Eq.  (4.10c)]. 

The  deflection  may  now  be  expressed  as 

2  G,  ... 

w=  [8(7A-?Ap)/e  J{  l-l  (7A-7Aj)/(M-7Ap)  ][git(x)  +  (g2(x)/4BT]]ri2)  (~  +  ^r)  Ja"1 

(4-15) 
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5-  STRESSES  IN  SHELL 


The  stress-displacement  relations  are  obtained  by  combining  Eqs. 

(A. 5),  (A -21)  and  (A. 23);  hence, 

ax  =  [Ex/(I-v2)][ex  -v  w/(a-z)  -  zuj  x!  (5 •  la) 

CT<p=  ly(l-v2)H-w/(a-z)  +  vx(p(eXo-2^x)]  (5-lb) 

From  Eqs.  (2.  I  la)  and  (2.17) 

exo  =  (,/ax)[Nx  +  axcp(w/a}  +  bx4°,x]  {5*2) 

It  should  be  noted  that  the  entire  analysis  to  this  point  can  be 
applied  to  a  shell  which  is  inhomogeneous  through  the  thickness-  In  an  actual 
filament  wound  composite  shell,  the  elastic  constants  are  discontinuous  functions 
of  the  thickness  coordinate  z ,'  hence,  the  actual  stresses  are  discontinuous 
functions  of  this  coordinate.  In  the  present  analysis  the  less  complex  state 
of  stress  in  an  equivalent  homogeneous  material  will  be  computed.  The  stresses 
in  the  actual  (nonhomogeneous)  shell  are  approximated  in  Appendix  E,  in  which 
suitable  multipliers  (stress  conce-n  t  rat  i  on  factors)  are  utilized. 

Equations  (5-1)  and  (5-2)  combine  to  yield,  for  a  homogeneous  material. 


=  1  - 


— (^w  -  (f+ff-M 

l-26(~)  h  5  h 


h  1+26  (-*-r)  (h 

f£l-  „  (|  (h/a)  r  1-v4  n 

Nx  V  V  1-26  $) 


-  ♦  irw 


where 


W(x)  =  E  w(x)/N 

A  X 

A  (x)  =  lExh2/2(l-v2)NxJu)(x)^ 


(5  -  3a) 


(5- 3b) 


(5- 3c) 


(5  3d) 


2,V 


The  shea'  stress  is  assumed  to  be  parabo  1  ica 1  I y  distributed  across  the 
shcl )  thickness.,  i  e. , 


7  =  J  (Q/h)  [  ? 

XZ  2 


(5-4) 


Thus,  the  maximum  shear  stress  occurs  at  the  median  surface.  At  the 
frame,  this  stress  is 


Txz  h  6g  A* 

"n  =  “lbgT  1  ^/?^A  "  A1J 

x  ^ 


(55) 


The  direct  stresses  can  be  separated  info  their  memb-ane  and  bending  components,- 
thus, 


membrane 


6A/5 


(5-6a) 


membrane 


(h/aJW/v^  -  5A/3] 


VN  ’ 

X 


(4?  V^W  +  A)  (2z/h) 


<p 


oh 

T2-) 

X 


bending 


t(a}  2z 

v[“frrw  +  A1(^) 

vv } 


(5 -6b) 


(5-  6c) 


(5-6d) 


in  which  the  approximation  that 


24. 


has  been  used.  A  shell  subjected  to  the  symmetric  type  of  load-support  system 
described  here  will  usually  develop  its  largest  stresses  either  at  irudbay  or 
at  the  frame.  Therefore,  the  computations  in  this  work  will  be  at  these  stations- 
Hence,  it  is  convenient  to  list  the  nond imens iona 1  parameters 


w  =  -  8(f)2  -4  [  — r - )  -  (z — r - )  a" J 

m  L  n  n 


=  A  -  K 


A  -  A, 


8(-)2(f)4e2E  §A_Ai 


~  /*  27  *  “  h  )(G4  -  2?G2)G  - 

(1-v  2 


f  h  A  -  An  £  A  -  A ,  ,v 

2  JL  r  (2 - )  .  (Z - )Alj 

*  U  h  '  '  h  ;G£j 

<P 


8(f)2(f)4e2E  ~  A  “  A  j  A* 

Af "  1  <c5  -  2?> t 1 

'■  v  ywcp 


(5 -7a) 

(5-7b) 

(5 -7c) 

(5  - 7d) 


where  the  subscripts  m  and  f  imply  that  the  function  is  calculated  at  the 
midbay  and  frame  shell  stations,  respectively.  The  stresses  can  be  calculated 
by  substituting  Eqs.  (5-7.)  into  (5-5)  [or  (5-6),  if  the  membrane  and  bending 
stresses  are  desired]. 

The  bear.-colurr-  term  y,  in  the  absence  of  transverse  shear  and 
initial  deflections,  has  been  shown  to  have  little  effect  on  the  stresses 
in  the  primary  region  of  interest  [see  Ref.  7J-  However,  the  combined 


effects  of  bcairrCuluftii  find  initial  dcjflcCtiGfi  tcfniS  &)  ,  Of  uf  beonr 

7 

column  and  transverse  shear  terms  (7^  8),  can  be  substantial.  This  can 
be  seen  from  an  examination  of  Eqs.  (.57)  and  (5-6).  For  example,  the 
bending  stresses  become  extremely  small  in  the  neighborhood  of  L  10  Aj  (they 


vanish  when  Flugge  terms  are  neglected).  When  transverse  shear  deformation 
is  considered  the  bending  stresses  may  be  considerably  affected,  depending 
upon  the  magnitude  of  ?• 
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6.  DESIGN  FORMULAE 


The  shell  stresses  and  deflections  may  now  be  utilized  to  develop 
convenient  formulae  for  several  quantities  which  may  be  of  interest  to  the 
des  i  gners  . 

A.  Circumferential  Normal  Stress  in  Ring 

Combining  Eqs.  (A. 26),  (3-2),  (3-19)  and  (5 -7c)  yield  the  follow¬ 
ing  expression  for  the  circumferential  normal  stress  in  a  ring: 


c  = 
s 


(E(pf/Ex)(Nx/h)lWf/(R  -  Z)} 
N  E_ 


32(f)2 

=  (4)  (-4)-—- 

e¥  <H-Z) 


2  A  -  A 


u2- 


-)  -  4 


i  A  “  A, 


-)a"j 


(6.1) 


B.  Equivalent  Line  Load  on  Ring 

If  a  line  load  Q.^  acting  at  the  ring  datum  circle  is  made  stati¬ 
cally  equivalent  to  the  total  inward  radial  load  acting  on  a  ring  segment, 
then  radial  equilibrium  gives 


Q. 


f 


-  F  /R 
s 


Then  Eqs.  (6.2),  (3.3),  (3-^)  and  (6.1)  combine  to  yield 

N  E  f  A,c'  „  -  A  "  Ao  -  A  "  A1 

Of  =  (i>(#)(-V)32(f) -  (-^ — K — )A“1 


where  Q.^. 


is  chosen  positive  in  the  inward  radial  direction. 


(6.2) 


(6.3) 


27- 


C-  Interaction  Pressure 

The  interaction  pressure  Rifiy  be  Calculated  f ron  E^s •  (5-20)  and 
(4. 15) •  Hence, 


P(l+6)2(ac/3*)[^f( 


l-Pc  A*  6«  _ 

j|)[(r/7)4  -  a2) 


(6.4) 


where 

A2  =  A0  “  L23cB[6*  -  (v^/2)(l+6)]/4a(l+6)(l-8c)G, 

-  \V/4a(l+6)3Hl  +  BBc/(l-ec)G,l  (6-5) 

Thus,  for  a  large  enough  outward  initial  deflection  (a2),  the  interaction 
pressure  vanishes  and  the  frame  would  be  ineffective. 

D-  Equivalent  Line  Load  on  Effective  Ring 

The  line  load  acting  on  the  section  of  shell  in  contact  with  the 
frame  is  defined  to  act  at  R  and  to  be  statically  equivalent  to  the  radial 
loads  acting  on  this  shell  section. 

%  =  (abc/R)(P6*  -  Pc(26“  -  (1+6)]  -  (2Qf/bc))  (6.6) 

where  pc  and  Q,^.  are  obtained  from  Eqs-  (6.4)  and  (4.15),  respectively. 

If  the  ring  and  the  section  of  shell  attached  to  it  are  considered 
together  as  an  "effective  ring",  the  line  load  acting  at  R  which  is  stati¬ 
cally  equivalent  to  the  applied  radial  loads  is  obtained  from 
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%**(  +  %  (6» 

E-  Axial  Location  of  Vanishing  Curvature  Change 

It  may  be  of  interest  to  the  designer  to  know  positions  along 
shell  generators  at  which  the  bending  stresses  are  negligible.  An  inves¬ 
tigation  of  Eqs ■  (5-6c)  and  (5.6b)  reveals  that  the  bending  stresses  be¬ 
come  higher  order  when  the  change  in  curvature  vanishes,  i.e.,  when 

A  =  0  (6.8a) 

or 

05  =0  (6.8b) 

i x 

Equation  (6.8b)  and  Eq  ■  (4.14)  then  combine  to  give  the  equation  to  deter¬ 
mine  the  point  x'  at  which  the  change  in  curvature  vanishes-  Hence, 


-(1/4Bti,ti2){[*(G4/-2)  +  1 3/2  +  S[(G4/G2)  +  *])g2(x') 

+  [(G4/2G2)  -  ?)g4(x')l  =  G  (6.9) 

x'  must,  of  course,  be  obtained  by  a  numerical  procedure.  It  is  noted  that 
when  beam-column  effects  are  considered,  x'  varies  somewhat  with  y.  There¬ 
fore  the  location  at  which  bending  is  minimized  (in  the  above  sense)  varies 
slightly  with  applied  pressure. 
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F ■  Effect i ve  Wi dth 

If,  as  in  some  ring  stability  analyses,  the  entire  structure  were 
thought  of  as  primarily  rings  with  modified  flanges  (i-e-,  the  shell  skin 
is  treated  as  an  additional  flange  on  the  ring),  then  it  is  of  interest 
to  establish  that  axial  length  of  shell  which  may  be  used  as  an  effective  ring 
flange  width-  Ihe  criterion  used  to  determine  this  "effective-width"  is 
(see  Ref.  5) 


1  1  1-/2  ,  , 

LN  (?)  =  f  N  (x) dx  =  IN 
e  qp'2'  .  <p'  q 


(6. 10a) 


average 


where 


Vx)  =  h<Wx) 


(6. 10b) 


L  is  the  effective  width-  is  the  resultant  of  the  modified  circum- 

e  <p 

ferential  membrane  stress  a_^(x)  which,  in  turn,  is  that  portion  of  the 
circumferential  membrane  stress  which  is  due  to  axial  bending.  0  *(x)  is 

<pr» 

obtained  by  omitting  terms  which  did  not  arise  from  axial  bending  considera’ 
tions  from  the  stress  given  in  Eq -  (5.6b).  Thus,  from  the  circumferential 
membrane  stress  one  must  subtract  the  hoop  stress  of  an  unsupported  cylin¬ 
der  under  hydrostatic  pressure  p;  this  leaves  as  that  part  of  the 

shell  membrane  stress  which  is  caused  by  the  ring  (which  causes  the  axial 
bending).  Hence,  from  Eqs.  (5-66),  (5-3c)  and  (4.15) 


<U(x>h  h  v  8E  7A  . 

- =  (j>  (r3*)  [w(x)  +  -f-2]  -  4  AW 

Nx  a  Vxtp  e2N  5 

T  X 


(O.11) 


30. 


Considerations  of  overall  equilibrium  of  a  shell  loaded  only  by  the  inter¬ 
action  pressures  leads  to  . 


L  N 

e  cp 


average 


=  Pcbc(l 


6  -  26") 


(6.12) 


Therefore, 


L  = 
e 


bcbch/Hx)  0“6"2$W)/[  ^vVh/vxq)a^  (Wf+8Ex7Ao/Nx02)-6Af/3] 


(6. 1 3a) 


1  I 

When  transverse  shear,  initial  deflection  and  Flugge  effects  are  neglected 
U=a=6=A=  0),  this  relation  reduces  to  faef.  7) 

L£  =  LG,  +  bc  (6.13b) 

G.  Axj svmmetri c  Buckling  Loads 

Axi symmetric  buckling  of  a  long,  unstiffened,  orthotropic,  circular 
cylindrical  shell  would  occur  either  at 


*  =  *cr  =  ' 

or 

B  =  B  =  0 

or 


(6.14a) 


(6. 14b) 


depending  on  which 
Appendix  C).  When 


of  the  above  conditions  occurs  at  the  lowest  j  (see 
:he  transverse  shear  parameter  §  is  relatively  small 


(G  large)  the  criterion  for  buckling  of  art  unsupported,  infinite  length 

X  2 


shell  Is  given  by  Eq.  (6.)4a)*  From  this,  the  axisymmetric  buckling  load 


i  s 


I  *  Vf 

- -  8  J  hE  D 

p’  =  (  . I  1  +  4|a  -  5  -  a)“5 - (when  1  -  4;r£  >0)  (6.15) 

*  a2 (1+6) 


The  ring-reinforced  shell,  due  to  its  greater  rigidity, 
possesses  an  axisymmetric  critical  load  p^  in  excess  of  that  of  an 
unstiffened  shell,  i.e.. 


p  >  p  or  7  >  7  or  4  >  1 

rcr  r  cr  '  ■ cr 


(6.16) 


At  the  buckling  load,  the  deflection  becomes  infinitely  large. 
Examination  of  Eqs.  (4.13)  and  (4.10b)  shows  that  this  occurs  when 


G  ( 1 
1  ' 


fi  1  [  G  " 
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=  n 
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However,  since  i)icr  >  1,  G|  in  its  present  form  [see  Eqs.  (4. lid),  (4.11c) 
and  (4-7) J  is  a  function  of  imaginary  quantities.  In  order  to  avoid  this, 
the  buckling  criterion  is  rewritten  in  the  form 


BP*  +  G5(1-PcMP*  +  0+6  -  26*)al/Pc  =  0 


(6.18a) 


where 


(6. 1 8b) 


is  obtained  by  setting 


T\\  =  -  iH,  =  .  *  “  l/(2/B) 


(6. 18c) 


in  the  equation  for  g  ^  . 

I f j  as  occurs  often  in  practice, 


P  <  <  ft  <  <  1 
rc  c 


(6. 19a) 


The  buckling  load  can  be  approximated  from  a  numerical  solution  for  i  of 


ac(26 


(6. 19b) 


If  more  accuracy  is  desired,  the  solution  of  Eq .  (6.19b)  may  be 
used  as  a  first  approximation  in  an  iterative  solution  of  Eq  •  (6.18a). 

For  a  rigid  frame  (shell  has  clamped  ends)  ct  -♦  <*>  .  Therefore, 
in  order  to  fulfill  the  buckling  criterion  (6.18a) 


G 


5 


—  CO 


(6.20a) 
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This  requires  that  [see  Eq .  (6.l8b)j 

B  i  i  B 

~V  sin  (T|.0)cos(T|,0)  +  cos  (T1  e)sin(T),0)  =  0  (6.20b) 

n,  i  ii2  2  z 

It  is  noted  here  that  the  results  presented  in  this  section  are  valid  for 
both  internal  and  external  rings. 


7.  YIELD  CRITERIA 


The  collapse  mechanism  of  filament  wound  composite  shells  is  not 
yet  fully  understood.  Experiments  indicate  that  collapse  occurs  near 
the  frame  before  the  buckling  pressure  is  reached  [Ref.  1  ].  Although 
the  severe  delaminat ion  noticed  in  ruptured  shells  indicates  an  ultimate 
shear  failure,  the  prevailing  opinion  seems  to  be  that  collapse  occurs 
after  a  complicated  sequence  of  the  combined  effects  of  yielding,  local 
instabilities  and  crack  propagation. 

It  is  not  the  purpose  of  this  work  tD  develop  a  mathematical 
means  of  describing  the  collapse  mechanism,  but  rather  to  describe  the 
stress  state  while  the  shell  material  behaves  as  an  elastic  continuum. 

It  would,  therefore,  be  advantageous  to  have  some  idea  of  the  range  of 
pressure  for  which  the  linear  elastic  laws  are  applicable- 

In  Ref-  II  ,  R-  Hill  developed  a  three-dimensional  yield 
criterion  to  be  used  on  orthotropic  materials-  He  introduced  a  "plastic 
potential"  function  of  the  stresses  which,  for  an  orthotropic  material 
in  a  state  of  modified  plane  stress,  with  transverse  shear,  can  be  writ¬ 
ten  in  the  form 


(7  1) 


where  f  is  the  "plastic  potential"  function,  Y  ,  Y  .  Y  and  Y  are  the 
r  r  x  cp  z  xz 

respective  axial,  circumferential,  radial  and  transverse  shear  yield 
stresses.  In  Ref.  7  >  in  which  a  composite  shell  theory  was  developed 
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without  introducing  effects  of  transverse  shear  deformation.-  the  correspond¬ 
ing  form  of  Eq.  (7-i),  i-e.,  t2X  =  0,  was  employed. 

The  yield  criterion  is  defined  to  be 


(7-2) 


Equations  (7-2)..  (7-0  and  (2-5)  combine  to  give  a  formula  for 
the  pressure  necessary  for  yielding  to  occur  at  a  particular  position  of 
the  shell.  This  "yield-pressure"  may  be  calculated  from 


Pv(*,z)  2(--)  o  h  o  h  V  o  h  V  „  o  h 

Yz  (l+6)2  Nx  NxnV  Nx  V  V' 


a  h  o  h  V  ,  t  h  ,  ,  /0 

1  KU  J  'Y  '  '  N  ' 

XX  XZ  X 


(7-5) 


When  beam-co 1 umn  effects  are  considered,  the  right-hand  side  of 
Eq.  (7-3)  is  dependent  on  p^  through  the  parameter  y.  In  such  a  case  the 
equation  is  transcendental  in  p^.  A  good  first  approximation  for  p^  may 
t;  obtained  by  taking  j  -  0  in  the  right-hand  side  of  Eq.  (7-3);  then,  if 
more  accuracy  is  desired,  the  yield  pressure  thus  obtaned  couid  be  utilized 
to  initiate  an  iterative  numerical  procedure  for  a  more  accurate  solution- 
Fiberglass  shel’s  are  generally  constructed  so  that  the  radial 
direction  is  normal  to  every  filament.  The  radial  yield  stress  (Y7)  of 
such  shells  is  the  relatively  small  yield  stress  of  the  resin.  Thus,  a 
good  approximation  to  the  yied  pressure  is  given  by 
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and 

{yy2  <  <  i  (7  -4c) 

The  most  severe  stress  situations  usually  occur  either  at  mid¬ 
bay  on  the  outer  shell  surface,  (x,z)  =  (0,  -h/2),  or  at  the  frame  on 
the  inner  shell  surface,  (x,z)  =  (±  L/2,  h/2).  The  yield  pressures  cor¬ 
responding  to  these  positions  may  be  calculated  by  using  Eqs .  (5-7),  (5-4) 
and  (5 -3)  in  Eq.  (7 -4a)  or  in  Eq.  (7  3)- 

An  estimate  of  the  pressure  P  (x)  at  which  yielding  has  pene- 

^o 

trated  through  half  the  shell  thickness  at  midbay  [see,  for  example.  Ref.  JL 
could  be  obtained  by  neglecting  the  bending  stresses,  i.e.,  by  using  Eqs. 

(3-4),  (5.6a),  (5-6c)  and  (5-7)  in  Eq.  (73)  or  (7-4) ■ 


The  manipulations  described  above  were  carried  out  in  Ref.  7, 
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to  result  in  the  following  yield  pressures  when  transverse  shear,  Flugge, 
and  initial  deflection  effects  are  neglected: 
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M  *  '  ‘  *  le'  +  -  fp*  j  <V| Cj  ♦  e5)j- 


py(o,  “  5)  =  3  to,  +  /ft  (7^)2  -  e.E.G.s.  +  0  F,2Gl  2i  .  A 


(7- 5a) 


V  ^.-2-4  +  Vi  VJ 


(7-5b) 


pyo(0)  =  (a)[^i  +  V«32a  +  ~  e5)]’,/2 


(7-  5c) 


in  which 


I ,  ,  6c6"A 

W('-ec)o,  ‘«c  •  rTTT^-1 


(7-5d) 


£-  =  2(^“^+7t) 

9  X  2 


(7-5e) 
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e  ,  =  - a— 
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Note:  for  approximation  (7-4) 


S  =  e5  =  2?1 


e  =  Vj£1 

C2  2 

2 


2  +  v. 


©4  = 


(7  •  6a) 

(7.6b) 

(7  ■  6c) 


The  simplest  failure  criterion  is  the  assumption  that  the 
material  fails  when  the  maximum  stress  reaches  the  uniaxial  yield  stress 
of  the  material.  Reference  7  also  applied  this  criterion  and 
arrived  at  the  following  results-- 
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8.  NUMERICAL  EXAMPLE  AND  DISCUSSION  OF  RESULTS 


Numerical  computations  were  performed  on  the  IBM  7040  computer 
located  at  the  Polytechnic  Institute  of  Brooklyn.  These  calculations  were 
intended  to  examine  the  separate  and  combined  effects  of  large  rotations 
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(beam-column  effect),  thick  shell  terms  consistent  with  Flugde  accuracy, 
transverse  shear  deformation  and  initial  deflections.  The  particular  glass- 
reinforced-plast ic,  ring-stiffened,  circular,  cylindrical  shell  of  interest 
was  one  previously  investigated  by  the  David  Taylor  Model  Basin  [Ref.  1 5  ]  - 
The  properties  of  this  structure  are  as  follows: 


The  shear  modulus  was  not  available  and  was  taken  to  correspond  to  similar 
structures  listed  in  Ref.  11. 


Geometry  (see  Fig.  l) 


a  =  3- 134" 


(continued  on  next  page) 


4!. 


The  results  of  the  above  calculations  are  presented  in  Tables  1 
and  2  and  Figs.  3  through  10.  The  calculations  in  (2)  and  (3)  were  prompted 
by  t.he  similarity  which  exists  between  the  isotropic  equations  found,  for 
example,  in  Ref.  (5)  (using  E  =  |  E}  and  the  orthotropic  equations 
[see  Ref.  7]-  In  the  few  cases  considered,  this  procedure  was  justified 
since,  for  both  the  thick  (h  =  0.388")  and  the  thin  (h  =  0.194")  shells, 
the  values  of  the  stresses  obtained  for  the  "equivalent  isotropic"  and 
the  orthotropic  cases  agree  to  within  10  percent  [see  Figs-  3  to  6  ]. 

The  stiff  ring  [E^-/E  =  1-335  in  Figs.  3  to  6  ]  caused  smal¬ 
ler  circumferential  stresses  but  higher  axial  stresses.  The  bending 
stresses  were  higher  for  the  stiff  ring  as  expected  [see  Figs-  7  -  10. ], 
and  the  corresponding  yield  pressures  could  be  slightly  lower. 

It  is  evidenc  from  Figs.  3  through  10  that  in  the  loading 
region  0  <  7  <  0-4,  the  nonlinear  beam  column  eFfect  is  small.  However, 
if  a  line  is  drawn  at  the  origin  tangent  to  the  curves  of  Figs.  3  through 
(10)  (yielding  the  linear  solution  which  neglects  the  beam-column  effect), 
it  can  be  seen  that  some  non  I ineari ties  first  appear  at  7  =  0-25  [e-g., 
see  Fig.  7  ]•  If  7  were  taken  larger  than  0.4,  the  non  1  inear i t ies  would 

no  longer  be  minor.  Of  course,  for  the  cases  considered,  when  7  >  0.4, 
p  >  37,000  psi  (for  h  =  0-388  in.)  or  p  >  10,000  psi  (for  h  =  0-194  in.) 
and  these  pressures  are  well  beyond  t.he  yield  pressures  that  were  obtained. 

The  yield  pressures  presented  in  Table  2  are  those  obtained  from 
Eqs.  (7-5),  (7-6)  and  (7-7)-  When  the  yield  pressures  were  calculated 
from  the  approximation  described  in  Eqs.  (7-4)  and  (7-6),  they  differed 
from  those  determined  by  Eq .  (7-5)  by  less  than  two  percent.  Thus,  the 
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low  yield  pressures  obtained  are  caused  by  the  smaii  radial  yield  stress. 

Hence,  if  the  resin  yield  stress  could  be  increased,  there  would  be  a  sub~ 
stantial  strengthening  effect  in  that  significant  increases  in  initial 
yield  pressure  could  be  realized. 

The  remaining  numerical  calculations  [Figs.  11  through  14]  are 
directed  at  investigating  the  effects  of  transverse  shear  deformation, 
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initial  deflections  and  thick  shell  (Flugge)  terms.  The  structure  con- 
sidered  is  basically  that  described  in  Eqs.  (8.1)  and  (8.2).  The  ring 
deflection  is  described  by  Eq.  (3-11). 

In  the  stress  vs  initial  deformation  curves  [Fig.  1),  13  and  17] 
the  bending  stresses  can  be  obtained  by  measuring  the  vertical  distance 
between  the  outer  and  inner  surface  stresses,  whereas  the  membrane  stresses 
are  equal  to  half  the  sum  of  the  outer  and  inner  surface  stresses- 

When  initial  deflection  and  the  beam-column  effects  are  negligible, 
the  membrane  stresses  exhibit  very  little  change  with  increased  transverse 
shear  deformation  effects  [see  Figs,  15  and  )6].  However,  the  axial  bending 
stresses  increase  at  the  frame,  and  decrease  at  midbay  The  lowest  (critical) 
yield  pressure  decreases  slightly  (See  Table  3)-  Although  in  general  the 
beam-column  effect  becomes  more  significant  when  transverse  shear  deformations 
are  considered  (see  Table  3),  this  effect  remains  negligible  for  pressures 
below  the  initial  yield  pressure. 

The  influence  of  initial  deformations  stems  from  beam-column  effects. 
Because  of  the  very  small  deflections  for  A  =  0,  referred  to  earlier  (see  Table  1), 
the  beam-column  terms  were  insignificant.  However,  with  initial  deflections  as 
large  as  A  =  ±  G.2h  such  terms  become  important.  Additional  deflections  (w)  con¬ 
tribute  very  little  to  beam-column  effects- 
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Thus.-  the  results  shown  in  Figs.  11  through  18  have  been  obtained  by  using 

(7/7)A  =  4  7=0  (8.3) 

in  Eqs.  (5-7) •  Numerical  results  for  7  ^  0  have  been  obtained  and  justify 
the  above  procedure.  For  example,  at  7  =  0.3  the  above  procedure  gives  the 
stresses  to  within  six  percent  of  the  more  accurate  value-  Hence,  instead 
of  plotting  numerous  identical  sets  of  curves  for  various  values  of  7,  it  is 
sufficient  to  present  the  single  set  of  results  obtained  by  using  Eqs.  (8.3). 

In  the  absence  of  transverse  shear  deformations  (|  =  0),  the  bend¬ 
ing  stresses  decrease  with  increasing  outward  initial  deflection  until,  at. 

4  =  A|  *  -  00/2  in.  for  the  shell  considered  (h  =  0.388  in.),  the  bending 
stresses  become  negligible  (Fig.  II)  and  the  shell  (with  A  =  A|)  behaves 
almost  as  a  membrane.  This  effect,  coupled  with  the  sharp  decrease  in 
hoop  stress  that  accompanies  large  outward  (negative)  initial  deformations 
causes  the  yield  pressure  to  increase  as  A  becomes  more  negative  (Fig.  12). 
For  A  <  Aj  yielding  can  begin  on  the  outer  surface  of  the  shell  at  the 
frame,  whereas  for  A  >  A]  yielding  begins  at  the  inner  shell  surface  of  the 
frame.  The  actual  initial  deflection  at  which  yielding  occurs  throughout  the 
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shell  wall  at  the  frame  varies  slightly  from  the  above  value  when  Flugge 
accuracy  is  maintained  (compare  Figs-  14  and  12). 

I  ( 

The  most  important  Flugge  type  correction  term  is  that  which  in¬ 
creases  the  axial  stress  resultant  Isee  Eq.  (2.5)]-  Further  comparison 
of  Figs.  11  and  13  reveals  that  the  stresses  (nondimensional ized  with  respect 
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to  N  /h)  are  only  slightly  influenced  by  Flugge  type  terms. 


The  actual  stresses 
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however,  are  proportional  to  and  change  by  about  10  percent  (for  h  =  0.588") 
due  to  the  change  This  shows  up  as  a  10  percent  decrease  in  yield  pres- 
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sure  when  Flugge  zsr  ,s  is  included  (see  Figs.  12  and  14) . 

The  effect  of  transverse  shear  deformation  is  to  increase  the 
bending  stresses  at  the  frame  and  decrease  those  at  midbay  (Fig.  15).  Also, 
the  value  of  the  initial  yield  pressure  (e  g.,  at  the  shell  inner  surface 
near  the  frame)  decreases  when  such  effects  arc  included  (Fig  16) ■  Comparison 
of  Figs.  11  and  17  reveals  that  the  changes  in  bending  stresses  described 
above  become  more  pronounced  with  increasing  initial  deflections-  Accordingly, 
the  initial  yield  p  assure  decreases  with  initial  inward  deflections  (Fig  18). 

Thus,  from  the  results  obtained  the  following  conclusions  may  be 

presented : 

1.  For  the  shell  material  considered  here,  an  equivalent  isotropic 
analysis  gives  stresses  accurate  to  within  ten  percent. 

2.  Stiffening  the  ring  may  increase  the  axial  bending  stresses 
enough  to  weaken  the  structure  by  causing  early  yield.  This  was  observed 
experimentally  in  Ref.  i. 

5-  The  relatively  short,  thick  shell  deflects  very  little  before 
yielding  and  is  not  susceptible  to  beam-column  effects.  Initial  deflections, 
however,  may  be  large  enough  to  cause  significant  beam-column  effects.  Short's 
suggestion  (Ref.  9)  that  the  beam-column  effect  be  utilized  to  strengthen 
the  shell  by  building  in  an  outward  initial  deflection  that  minimizes  the 
bending  stresses  may  be  used  on  this  type  of  shell.  All  that  need  be  done 
is  to  wind  the  shell  on  a  barrel  shaped  mandrel. 
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k.  Yield  pressures  which  are  well  below  experimental  collapse 
pressures  indicate  that  yielding  begins  well  in  advance  of  collapse.  This 
may  be  remedied  by  increasing  the  yield  stress  in  the  radial  direction. 

5-  If  the  shell  is  permitted  an  additional  degree  of  freedom  by 
allowing  transverse  shear  deformation,  the  frame  bending  stresses  increase 
and  the  initial  yield  pressure  decreases.  The  transverse  shear  deformation 
theory  cannot,  however,  be  considered  as  necessarily  more  accurate;  this  is 
demonstrated  in  Ref.  2- 
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6.  The  only  Flugge  type  term  that  has  a  significant  effect  is 
that  due  to  using  the  outer  shell  diameter  when  calculating  the  axial  pres¬ 
sure  force  acting  on  the  shell. 

7.  For  the  shell  structures  examined  in  the  present  work,  the 
yield  condition  was  reached  before  y,  £,  q  or  6  exhibited  any  substantial 
influence  on  stress  curves.  Thus,  it  was  found  that  a  perturbation  technique 
can  be  employed  to  give  fast  and  accurate  results  (see  appendix  0) . 

8.  In  order  to  obtain  reasonable,  estimates  of  the  stresses  in 
the  constituent  materials  it  is  necessary  to  modify  the  stresses  obtained 
for  the  equivalent  homogeneous  shell.  This  modification  can  be  achieved 
by  multiplying  the  stresses  given  by  Eqs-  (5*3)  by  appropriate  stress  con¬ 
centration  factors  (see  Appendix  E) . 
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APPENDIX  A 

BASIC  SHELL  EQUATIONS 

The  differential  equations  describing  the  axisymmetric  deforma 
tion  of  the  shell  are  obtained  by  utilizing  approximate  forms  of 
the  basic  equations  of  the  theory  of  elasticity. 

First  the  strain  components  are  found  in  terms  of  the  displace 
rr.«.ts  and  rotation-  The  length  of  an  infinitesimal  line  dsQ  in  the  well 
of  a  circular  cylindrical  shell  with  a  slight  axi symmetr ic  initial  de¬ 
flection  w  may  be  obtained  from  (for  example,  see  Ref.  6) 
o 

dsQ2  =  d£o2  +  (a  -  z  -  wQ)2d<p2  +  d£Q2 
where 


*  x 


zw 


o.x 


r  =  z  +  w 

o 


w  =  w  (x) 
o  o'' 


(A.l) 


x,  9  and  z  are  Lagrangian  coordinates  in  the  axial,  circumferential  and  in¬ 
ward  radial  directions,  respectively,  and  aie  measured  from  a  point  On  the 
median  surface  at  the  midbay  station  of  the  shell  (see  Fig.  1).  These  co¬ 
ordinates  are  fixed  to  the  shell-  E  and  r  are  Eulerian  coordinates  in 

G  O 

the  axial  and  inward  radio)  directions,  respectively,  measured  from  the 
median  surface  at  the  midbay  station  of  a  perfect,  undeformed  shell;  a  is 
the  radius  of  the  median  surface  and  is  the  initial  inward  deflection  of 
that  surface.  A  comma  followed  by  a  subscripted  variable  denotes  differen¬ 
tiation  with  respect  to  that  variable- 


**9- 


After  ax i symmetric  deformation  the  infinitesimal  lire  attains  a 
new  length  d»  |  where 


ds  2  =  d£.2  +  (a  -  z  -  w  -  w)^dq/  +  d£.4 


2  .2 


I  M 


(A-i) 


wi  th 


S,  =  x  +  u 
C,  =  z  +  wQ  +  w 
u  =  u(x,z) 

w  =  w(x) 

§ ,  and  Cl  are  Eulerian  coordinates  similar  to  Z  and  C  .  u  is  the  axial 
1  1  o  too 

displacement  of  a  point  in  the  shell  wall  and  w  is  the  inward  deflection  of 
the  median  surface  from  the  initially  deformed  shape.  If  normals  to  the  median 
surface  remain  straight  and  unextended  throughout  the  deformation, 

u  =  u  -  z  (ui  +  )  (A  .  5a) 

Of  A 

u>  “  «  „  -  (A  Ob) 

;  ?■  A*. 

where  li  is  the  axial  displacement  of  a  point  in  the  median  surface,  u>  is  the 
additional  rotation  of  the  median  surface  normal  about  a  circumferential  axis 
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due  to  loading  and  y  is  the  corresponding  transverse  sheading  strain-  It 

is  noted  here  that  y  and  ui  are  not  independent  of  each  other. 

The  strain  components  are  obtained  from  [see  F.qs.  \A.l)  and  (A  T)] 


ds,2  -  dso2  =  2exdx2  +  2e^(a-w(j-z)2dcp2  +  2ezdz2  +  7yzciK<iz 
Hence, 


e  =  u  +  —  w  (w  +  2w  )  -  z  (w  -  y  ) 
^x  ,x  2  ,x  ,x  o,x  ,xx  xz,x 


e<p  11  "  w/(a-z"«D) 


e  =(w  “7  )w  ~  0 

z  ,x  x z  o,x 


(A.<0 


y  =  w  -on-uw  »  w  -  oo 
xz  ,x  ,x  o,x  ,x 


in  which  it  has  been  assumed  that 


(A.  5) 


<  <  w 


(w 


2w  ) 
o,x 


<  < 


w<<a  -  z  -  w_ 


(A.  6) 


It  is  noted  that  a  transverse  shear  term  has  been  retained  so  that 


i 


the  usual  Kirchhoff  assumption  that  lines  originally  normal  to  the  median 
surface  remain  normal  is  removed.  However,  with  regard  to  kinematic  con¬ 
siderations,  it  is  assumed  that  such  lines  remain  straight  and  uncxtended 


>1. 


In  order  to  obtain  a  set  of  equilibrium  equations  consistent  with 
the  strain-disp  lacement  relations  [Eqs-  (A. 5)]  the  frsl  variation  of  the 
total  potential  can  be  set  equal  to  zero 

6(11  +  V)  =  0  or  611  +  6V  =  0  (A. 7) 


where  U  is  the  strain  energy  and  V  is  the  potential  due  to  externally  apolied 
loads.  Thus, 

L/2  h/2 

6U  =  27r-l-L/2  ■('■h/2  (ox6ex  +  °<p%  +  Vez  +  Txz^xz)(a-z-Wo)dzdx  <A'8) 


When  Eqs.  (A. 4)  and  (A. 5)  are  substituted  into  Eq .  (A. 8)  and  integrations  by 
parts  are  performed,  Eq .  (A. 8 )  becomes 


L/2 


6U  =  27ra.f.L/2  U-NX;X  +  (Qwq.x)>x]6u  +  t-M^vx  -  <NX«  x^xw0, x>  "  f  ]*W 


L/2 


+  ("Mx,x‘Hl)5yx^dx  +  ]-L/2 


L/2  L/2 

+  [^x,x+  NxV)6wJ-L/2  +  [MxKz  “  6w,x  )1  -L/23 


(A.  9) 


in  which  the  stress  resultants  are  defined  by 
h/2 

N  =  f  o  [1  “  (2z/h) 6  -  (w  /a)]  dz 
x  v_h/2  X  o 

h/2 

N_  =  f  o  dz 

cp  J  .  ...  ffl 

-  .h/7  * 


h/2 

M  =  f*  o  [I  -  (2z/h)6  -  (w  /a)]zdz 
*  _u  /o  X  0 


h/2 

Q  =  K  r  T  l 1  -  (2z/h) 6  -  (w  /a)]  dz 
-h/2  ° 


(A.  10a) 


\t\  »  luoj 


(A.  10c) 


(A.  lOd) 
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where  0  is  the  axial  stress,  a  is  the  circumferential  stress  and  t  is 

x  CP  y.z 

the  transverse  shear  stress  with  k  as  the  corresponding  shear  coefficient. 
The  potential  due  to  external  pressuie  is  obtained  by  considering  the 
volume  change.  Thus 


(continued  on  next  page) 
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+  2C  +£ 


—  W 

W  ,  ,  X 

2a  a 


(6u 


+  j  6ui))dx 


+  p7raZ[  (1  +  27  “  “)2  (*u 


+  “  6o>)] 


L/2 

>1/2 


(A. 14) 


Finally,  the  variation  in  the  potential  due  to  other  externally  applied  loads 
is  given  by 


5V,  =  [ -N  5u  ] 


L/2 

-L/2 


-  [N  w  6w] 


L/2 

■L/2 


xz 


,)] 


L/2 

-L/2 


L/2 

”  [Mx,x6w]-L/2  +  lVo,x6u 


] 


L/2 

-L/2 


(A.  15) 


The  potential  of  externally  applied  loals  Is 


V  =  V  +  V, 
P  f 


(A. 16) 


Equation  (A. 7)  to  (A. l6)  leads  to  the  differential  equations  of  equilibrium 


-  N^x  +  p[i  +  h/2a  -  (vHnvo)/a](v^vo)  ^  >K  =  0 

”  ^x^o^x^x  ’  V3  '  Q,x  •  p[l  +  h/2a  "  <^o)/aJ(1  +  u,x 
M  -  Q  +  ^  [  !  *  ii/2a  -  (w  -r  w  )/aj(w  +  w  )  =  0 

j  JK  4  O  v  ^  X 


(A. 17a) 


+  h/2  a>  )  -  C 

(A- 17b) 
(A. 17c) 


with  boundary  conditions  at  x  =  ±  L/2 


N  -  Qw  +  (pa/2) l 1  +  h/2a  -  (w,+w  )/a]2  =  0  or  6u  =  C  (A. )8a) 

X  O;  X  I  O^- 

-  +  (n/2a)  (pa2/2) [ 1  +  h/2a  -  (w  +  wQ)/a]2  =  0  or  =  0  (A. l8b) 

N  (w  +  w  )  +  Q  =  C  or  &w  =  0  (A.  1 8c) 

x  o  ,x 

Equations  (A. 16)  combine  to  yield  the  resulting  axisymmetric  equilibrium  equa¬ 
tions  for  hydrostatic  pressure  loading  of  a  shell  of  uniform  thickness  (in 
which  the  assumptions  given  by  inequalities  (2.4)  have  been  utilized),  viz., 

-  N  +  2p[ 1  +  6  -  (w  +  w  )/a](w  +  w  )  =  0  (A. 19a) 

A  ^  A  U  U  y  A 

lNx(w  +  wo)^]<x  +  Mx>xx  +  N6/a  +  p[l  +  S  +  (w  +  wo)/a)  =  0  (A.  19b) 

-  M  +  Q  =  0  (A. 19c) 

X  f  X 

1 1 

where  6  =  h/2a  (set  equal  to  zero  when  Fluggc  accuracy  is  not  required,  ie-, 
when  the  variation  of  radius  of  curvature  across  the  thickness  is  neglected). 

For  the  problem  being  considered  here,  the  following  conditions  must 
also  be  satisfied  at  the  shell  boundaries  x  =  ±  L/2: 


N  +  (pa/2) [ ]  +  6  -  (w  +  w  ) / a ] 2  =  0 

A  U 


(A. 20a) 


w  =  w 

s 


(A. 20b) 


u>  =  0 


(A. 20c) 
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where  w„  is  the  deflection  of  the  ring  at  the  contact  surface  and  w  is 
R  s 

the  deflection  of  the  section  of  shell  in  contact  with  the  ring.  For  the 
ring-supported  shell  considered  here  the  axial  stress  resultant  at  the 
frame  is  due  to  the  externa!  pressure  loading  [Eq.  (A.2Qa)];  the  ring  and 
the  shell  remain  in  contact  throughout  the  deformation  [Eq.  (A. 20b)],'  and 
since  the  ring  flanges  are  not  permitted  to  bend, the  rotation  vanishes  at 
the  boundary  [Eq.  (A. 20c)]. 

The  stress-strain  relationships  are  those  for  an  orthotropic, 
I'-neariy  elastic  shell  with  se  i^ces  of  elastic  symmetry  defined  by  normals 
directed  along  the  coordinates  x,  <p  and  z,  respectively.  The  shell  skin  is 
assumed  to  be  in  a  state  of  modified  plane  stress  (i.e-,  with  a  transverse 
shear  term  permitted),  thus 


a 

x 


7$ (e*  *  W 


(A .  2  ia) 


o 


9 


u 


9 


Vxcp€x^ 


(A. 2  It;) 


T 


XZ 


G  7 
xz  xz 


(A- 21c) 


o  - 
z 


0 


(A. 2  Id) 


and  1  rom  I. be  axi symmetric  nature  of  the  structure  ar.d  loading 


T  =  T 

xq>  qtz 


0 


(A . 2  le) 


whc.  <*e 


v  E  =  v  £ 
x<j>  9  <jfc  x 


(A  ’If) 


V  V  =  V 

xcp  qx 


(A  •  2  ) g) 


i  n  wh  i  r  h  F  H  F  a  ro  t  lao  *»V  !  1  1  -ir\rl  /-*«•  —  .  ■  pm  1  _  1  -  _  »  ■  _J  »  ■ 

. .  - -  w.w  W..S,  «">  •  u  «  qny  VMWUIMICICIIUC!  Ct  i  d  :>  U  «  C  IHUUUIIj  rtib" 
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pectively,  G  is  the  transverse  shear  modulus,  V  and  v  are  the  Poisson 

xq>  qx 

rat i os ■ 

In  order  to  formulate  the  problem  in  terms  of  the  stress  resul¬ 
tants,  rather  than  the  stresses,  Eqs .  (A.10)  and  (A. 21)  are  combined.  Thus 


Hx  =  [ax  “  \(wo/a)]sx  +  Wq>  “  [bx-Bx(wo/a)]w, X  (A-22a) 

o  T  To 


=  t[  1  -  H/a)  M  -  [1  -  2  (w  /a)  J  B  +  Cl.  +  a  c  -bo, 
o  q)  o  9  9  x<p  x  j  j 


(A. 22b) 
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h  =  [  b  -  B  (w  /a)  ]€  +  fc 

x  x  x'  o  '  x 

o 


^  q> 


-  [dx  -  C  (w  /•)]«, 


x  o 


w 


Q.  =  C  (w  -  <ju)  -  (u)  -  x) 

Xx  .X  |  ,x 


in  which 


e  »  u  +  (w  /2) (w  +  2w  ) 
xQ  ,x  ,x  ,x  c,x' 


and 


h/2  E 

\  *  J . ..  r'1? 


-h/2  1-v 


h/2  E 

A  J  dz 

-h/2  l-v 


h/2  m 

b  =  J 

<P*  J  l -  .  2 


zdz 


-h/2  l-v 


J 


h/2  E  v 


*Jc* 


J -h/2  l-v2 


7T  dz 


h/2  E 

bx  =  J  “"*2  (l"  7)z?lz 
-h/2  l-v  3 


Bx  =  J' 


h/2  L 


x  J  u/,  ,  2 

-h/2  l-v 


zdz 


H/2  En-v 

d2 


,  h/2  E 

B  =  ~  f  zdz 

9  3  v -h/2  l-v2 


i  h/2  E  , 

c<p  tJ  -t  z  qz 

a  -h/2  l-v 


h/2 

r 

-h/2 


a  .  2 

I  -v 


C  =  f 

x  J 


XZ 


h/2 

m 

N 

r 

X 

-h/2 

l-v2 

h/2 

G  ( 

-h/a 

X  2 

dz 


(Av22c) 

(A . 22d) 

(A. 2» 

(A. 23b) 


(A.  24) 
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The  basic  Eqs.  (A.5)>  (A. 19)  and  (A. 22)  subject  to  boundary  con 

ditions  given  by  Eq.  (A. 20)  describe  the  deformation  of  the  shell.  These 

consist  of  ten  equations  in  the  unknowns  Q,,  yt,  u,  to,  e^, 

and  y 
xz 

Integration  of  Eq ■  (A. 19a)  and  application  of  the  boundary  con¬ 
dition  Eq.  (A.20e)  reveals  that  is  very  nearly  constant.  In  fact,  the 
assumption  that 


(w  +  wQ)/a  <  <  1  +  6 


(A. 25) 


leads  to 


Nx  =  -  (pa/2) (1  +  6)2  (A. 26) 


If  Eq.  (A. 22a)  is  used  to  eliminate  e  from  Eq.  (A.22b-d)  and 

o 

if  Eq.  (A.2Jb)  is  introduced,  the  resulting  equations  reduce  to  the  linear 
relations 


M  =  z  ,N  -  63E"w  -  0*0) 

x  lx  J  qix  x  , x 


N_  =  vlN  -  26E "w  -  53Elto 

H*  A  y  s  ^ 


Q.  *  hG  [  1  -  (w/a)  ]  (w  -o>) 

A*.  i  A 


=  hG  [I  -  (-8-)] 7 
xz  a  xz 


(A. 27a) 


(A- 27b) 


(A. 27c) 


where 


5i'- 


2,  =  [bx  -  Bx(wo/.)]/[ax  -  A  (w0/a)J 


6^E  =  b  -  a  z , 

5  q&  XCp  1 


0  =  d  -  C  —  -  [b  -  B  (w  /a) ]z , 

x  x  xa  x  x'ol 


v*  =  -  Vfr  - 
<p*  w 

.  o 
a  -  A  — 
x  x  a 


%  =  "  axV  +  1  1  “  (wo/,)]A9  "  M  “  2(wo/a)  %  +  C<p 


hG*  =  C 
xz  xz 


(A. 28) 


In  view  of  the  fact  that  Nx  Is  a  known  constant  given  by  Eq.  (A-26j) ,  the 
three  remaining  stress  resultants  are  seen  to  be  linear  functions  of  w 
and  o>- 

In  order  to  avoid  variable  coefficients  in  the  final  differential 
equations  the  shell  material  must  be  assumed  to  be  homogeneous  in  the  axial 
direction;  in  addition  the  w0  terms  must  be  neglected  in  Eqs.  (A-28)-  Thus, 
must  be  small  enough  for  the  following  assumptions  to  apply: 


w  /a  <  <  1/2 
o 


V<  <•*,  -  2w/(c*  -  zie*>  ■  'o 
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w 

-2-  <  < 


0(1) 


w 

-a 

a 


<  < 


0(1) 


(A.  29) 


If  the  shell  material  is  also  homogeneous  through  the  thickness,  the 
condition  wo/a  <<  1/2  prevails  and,  hence  Eqs .  (A. 28)  reduce  to 


*,  =  -  (a/J)62 

E*  =  v  E  /(l-v2) 

<px  x  ' 

°*  -  0X[ 1  -  (63/3) ] 
e*  =  V 1  +  (»3/5)/o-v2)'j 

Vcpx  ~  v<px 

*  i  i 

,  ,  t 

G*  =  kG 

XZ  XZ 

where 

0x  -  (l/!2)Exh3/(l  -  v2) 


(A.  28a) 
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When  £qs-  (A. 27),  (A. 28a)  and  (A. 29)  are  substituted  into  Eqs . 

(A. 19b)  and  (A. 19c)  the  shell  problem  reduces  to  the  solution  of  the  follow¬ 
ing  pair  of  linear,  second  order  differential  equations  with  constant  co¬ 
efficients: 


Lllw+  L12a3,x 


LI3Wo  +  L1 4 


(A- 30a) 


L2,w 


L22% 


=  L2JWo 


(A. 30b) 


where  L. .  are  cummutative  differential  operators  defined  by 


L  =  (-N  +  ~  ~  <-N  +  %  ~ 

11x3  cpx  4  '^2  a  cp  —  '  x  3  <pK  ^v2  a  9 


dx 


L12  =  D  ^2  +  ¥■  65E* 
12  x  dx2  ’ 


'21 


-  (Z*  *3C*  ^  ur.*  .  Ph>d  /2a  3C*  .  Kr*  ,d_i_ 

V  X  w  "»  MU  T  n  J  ft  ^  (  4  6  t  ^  hG  )  A 

3  9x  X2  2  dx  ^  ^  XZ  Jv2 


L 

dx" 


...  .2 

L,  ■  =  D  *-7  -  hG 
2c  x  .  I  xz 
dx 


L ,  1  »  (N  +  “r-  ~  N  — -7 

*  4  dx*  “  x  dx 
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Li4 =  ~r~Nx *  +  6> 


1  =  E& 

23  2  dx2  - 


;he  approximations  indicated  above  are  associated  with  terms  which 

I  l 

are  due  to  rruplecf  beam-column  and  Flugge  effects.  For  the  shell  considered 
here  these  terms  arc  extremely  small  and  hence  are  neglected. 
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Elimination  of  co  from  Eqs .  (A.JO)  yields  the  following  fourth  order 
equation  in  the  deflection  w: 

»,,„x  *  Me/i>V<»w)XX  ♦  Mo/iA/i  -  (M/ 8)%xxxx 

-  [47(e/L)2/B]w  +  8(e/L)2T[l  -  (v*  V2)(l+6)]/aB(l+6)  (A. 51) 

Q  j  X  X 

where 

A  =  y  -  I  +  a 

B  =  1  -  4  7i 
(9/L)U  =  hE*/4a2D* 

7  =  7  =  -  Nv/4oV(e/02 

X  X 

I  =  (0/L)2D'  /hG 

x  xz 

a  =  a63E*/5D*(0/L)2  (A. 52) 

x 

Subtraction  of  Eq .  (A - 50b)  from  Eq.  (A.30a)  yields  a  convenient 
expression  for  the  change  in  curvature 

«  =  [Bw  -  4|(0/L)2w  +  8*[1  -  (v”  /2)(l+S))/*(l+8)-47tw  l/(H2a6) 

y  X  f  XX  O  y  XX 


(A.  53) 
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Eqs.  (A. 19c),  (A- 27a)  and  (A. 33)  combined  to  give  the  following 
relation  for  the  transverse  shear  stress  resultant: 


Q  = 


D“{w 
x  ,xxx 


-  2(e/L)  (2f-ev)w 


4/fw  V(l+2or5) 
o,xxx  ' 


(A. 34) 


Finally,  Eqs.  (A. 34)  and  (A. 27c)  combine  [using  conditions  (A. 29)]  to  give 
the  rotation  as  a  function  of  w 


O)  =  {[1  +  4$(<y-01w  v  +  l('<-/e)2(Bw-4nw)  *  l/'.l+2o&)  (A. 33) 

The  governing  differential  equations  of  the  shell  consist  of 
any  two  of  Eq.  (A-3l),  (A  33)  and  (A. 35)- 

In  Eqs.  (A. 32)  the  nondimens ional  parameter  y  is  a  measure  of 
the  nonlinear  "beam-column"  effects  (see  Ref.  5)-  Although  y  and  7  are 
numerically  equal,  y  terms  result  from  nonlinear  effects  while  y  terms 
result  from  the  linear  effects  of  the  pressure.  Thus  if  the  beam-column 
effect  is  neglected  y  -  0  while  y  ¥  0. 

The  nond imens i ona I  parameter  £  is  a  measure  of  the  transverse 
shear  effect.  If  no  transverse  shear  is  permitted  G^z  and  |  ->  0- 
Furthermore,  the  nond imens i ona!  parameters  a  and  5  vanish  when  Flugge 
type  accuracy  is  not  desired. 

It  is  noted  that  special  cases  exist  when  A  or  6  vanish.  When 
5=0,  y  =  0,  A  =  0  and  B  ■  1,  Eq.  (A-31)  reduces  to  the  equation  correspond¬ 
ing  to  the  omission  of  both  transverse  shear  and  beam-column  effects-  When 
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the  transverse  shear  and  beam-column  effects  combine  to  cau.c  G  =  0,  the 
entire  nature  of  the  differential  equation  changes.  When  the  higher  order 
effects  are  omitted,  i.e.,  when  7=§  =  6=  CX=A=0,  except  the  6  term 
in  the  denominator,  Eg.  (A. 31)  corresponds  to  that  obtained  in  Ref.  1,  in 

I  I 

which  Horn  provided  for  the  only  significant  Flugge  type  term,  which  cor¬ 
responds  to  6  4  0  in  the  expression  for  the  axial  stress  resultant  [see 
Eq.  (A. 26)]. 

When  transverse  shear  deformations  are  neglected  (£  =  0)  only 

one  equa.ion  is  necessary  to  describe  the  shell  deformation.  In  this 

case  Eqs.  (A. 33)  and  (A. 35)  are  satisfied  by  the  condition  that  u)  =  w  ; 

> x 

and  Eq-  (A-30  reduces  to 


w  wvvw  +  Me/L)2(7+a)w  +4(e/L)\y 

; AaAA  j AA 


-  47(6/L)2w 


o,xx 


+  8(6/L)27l  Mv"  /2)(l+t)]/a(!+6)  (A. 3U) 

<px 

l  I 

If  Flugge  thickness  terms  are  also  neglected  (a=5=0)  and  if  no  initial 
deflections  are  considered  (wq  ?  0),  Eq.  (A.3la)  reduces  to 


w,  xxxx+i+7^e/  L)  2w,  xx-*^  (®/L)  ^8  (8/  L) 2  7  ( 1  -v^/2)  /a 


(A- 31b) 


Equation  (A-3lb)  corresponds  to  that  obtained  in  Ref.  J,  in  which  similar 
formulae  as  those  presented  here  w  re  developed  and,  in  addition,  the  re¬ 
sults  were  expressed  explicitly  in  terms  of  the  load  and  the  basic  shell 
material  and  geometrical  input  quantities- 
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APPEND iX  6  (Ref.  7  ) 
EVALUATION  OF  c' 


The  geometric  ring  parameter  c’  is  defined  as  [sec  £q.  (3.4)] 


h. 


±  f 

A,  J 


Lliil 

h  R  -  Z 


dZ 


-S.  p 

Af\. 


C 


o  i 

where  Q  =  R  -  Z.  From  the  definition  (Bl),  it  can  be  seen  that  c* 
is  dependent  only  upon  the  geometry  of  the  ring  cross  section  and  is 
therefore  independent  of  the  choice  of  origin  of  Z- 

Performing  the  integration  indicated  in  (Bl)  for  a  general 
section  yields 


(Bl) 


"T 


I  I 


c' 


—  t.  t  . 

jf-  £(b.-hw)log(l  +  £-4  -  (bo-hw)log(l  -  -  hJogO  -  £-))  GENERAL  <V' 

f  i  co 

(B2) 


where,  d  =  Rq  -  R.  is  the  total  height  of  the  ring  cross  sectior  • 

A-  Origin  of  Z  at  Mid-Height  of  Rina: 

For  this  choice  of  the  origin  of  Z 


Therefore,  Eq.  (B2)  becomes 
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?  ti/R 

s  T  i(bj-hw>,03M  +  — "~I 

Af  '  w  1  ~d/2R 


<bQ-t,Jiog[  i- 


°/R  1  -  l-<*/2R  n 

- -  h  i 031  —  r:J] 

l+d/2R  l+d/2R 


GENERAL  "i' 


(83) 


However,  recalling  that 


a.'  .  II 

I og ( 1  +  x)  =  E  (“l)n  “  for  -)  <  x  <  1 

n=l 


log(l-x)  = 


a>  n 

2  s- 

n=  1  n 


for  -1  <  x  <  I 


and  that 


e 

n^O 


(-x)r 


for  -1  <  x  <  I 


permi ts  Eq-  (B3)  to  be  rewritten  as 


i  ^  R  r  /.  .  s  i  r  i  .  1  / 1  I  1  ^  1  .  d  £  _  U  i 

c  ~  r~  i  (8:-h)—  U  +  2  (•  2  ^ 

Mf  ,WR  £RR  R  R'R 


(b0-hw)^M  -}(' -J=)"+ !=?■(' 

~  RR  R 


R  ''  R 


+  hw  =  f  *js  <=>2)) 


GENERAL  "I" 


(B4) 


in  which  all  terms  of  third  or  higher  order  have  been  neglected  due  to  the 
fact  that  t  and  d  are  small  compared  to  R,  ie-. 


i  -  « i  -  (i)3  ~  i 


where  here,  and  in  subsequent  expansions,  R  is  considered  tc  be  very 


much  greater  than  d. 


I-  Rectangular  Section 


For  a  rectangular  ring  cross  section  of  height  d  and  width  h  , 


A,  =  dh  and  t .  =  t 

f  w  10 


Therefore,  Eq -  (B4)  becomes  (R  >  >  d) 


C '  rect 


RECTANGULAR 


2-  Equal  Flange  'T"  Section 


For  an  equal  flange  ‘T’1  ring  cross  section, 


b  =  b.  ='  b 
o  i 


t  =  t.  =  t 
O  I 


Therefore,  Eq .  (B4)  becomes 


‘•i  -  a7  *  12  *  f  ’  * *  *  f'- +  1  *  ilr 


ECtUAL  "I" 


where, 
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A£  =  2l(b-h  )  +  dh 
f  w  w 

If  the  second  order  terms  in  t  and  d  of  (b6)  are  dropped  as  being  small 

(R  >  >  d) 


For  an  idealized  "1"  section 


Therefore,  Eq .  (b6)  becomes  (R  >  >  d) 

c!  *  1  +  r  (-)2  ~  1  IDEAL  ’1" 

i  4  n 


A,  A,  Z  dAf 

;i  • 1  *7 1 2dAf  +  rri 

ft  R  1-Z/R 


However,  since  the  origin  of  Z  ia  taken  to  be  the  centroid  of  the  cross1 
sectional  area. 


Therefore, 


where. 


ZdAf  s  0 


c*  =  1  +  (-)  c 
ft 


p  =  V  I  f/A 


is  the  radius  of  gyration  of  the  cross-sectional  area, 


2 

1 =  f*  Z  dAf  is  the  moment  of  inertia  of  the  cross- 


sectional  area  about  Z  =  0, 


=  _L  r  Lzf  =  _L  r 
Tf  *  i -z/pT  If  - 


Z  b(Z)dZ 

h  1  -  Z/R 
o 


c 


70. 


Expanding  Z  /(I  -  Z/R)  in  a  series  (R  >  >  d),  performing  the 

integration  indicated  in  (B9)  for  a  genera)  "Il!  section,  and  neglecting 
5  —5 

terms  of  order  (Z  / 9r)  and  higher,  yields  the  following  expression  for  c 


_L  r  bo 


8  '  ^  V-.) 


(h  -t  )4]  +  l  (h.-t.,  -  {h  -t  )Uj)  GENERAL  "I"  (B 10) 

4r 


7  !• 


APPENDIX  C.  AX  I  SYMMETRIC  BUCKLING  LOAD  OF  INFINITE,  UNSUPPORTED 
ORTHOTROPIC  CYLINDER  UNDER  AXIAL 
COMPRESSION 


The  differentia)  equation  governing  the  deflection  of  a  perfect 
shell  subjected  to  axia!  compression  is  taken  from  Eq.  (2.1);  i.e.. 


M^)2  (7-5*3) 
w,  xxxx  +  1  -  4/5  xx  +  1  -4  75  w  -  0 


(ci) 


The  nontrivial  (ax  i  s  ymine  t  r  i  c  buckling)  solution  is 


w  =  DcosX 


(C2) 


where  L  is  the  wavelength  of  the  buckling  pattern.  Substitution  gives 


4  2  2  k 

X  -  40  ax  +  40  b  =  0 


(C3a) 


1-475 


b  = 


1-471 


(C3b) 

(C3c) 


Thus 


X2  =  202  (a  ±  a2-b) 


(C4) 


In  order  for  w  to  remain  finite  X  must  be  real-  Then  the  following  two 


poss ibi 1 i t i es  exist: 
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0  <  b  <  a^  (C51*) 

or 

b  <  0  (C5b) 


In  the  first  case  substituting  Eqs.  (Cjb)  and  (C3c)  into(C5)  yields  an 
expression  for  the  buckling  load  parameter  y 


/'  =  { 


1  +  kia  -  I  - 


O 


(C6a) 


when 


(c6b) 


For  sufficiently  small  a  Eq.  (C6a)  is  closely  approximated  by 


7  =1-0-1 (l-2o) 


(C7) 


In  the  second  case  (b  <  0)  the  buckling  load  given  by  Eq .  (C^b)  is 


7\ 


(C8) 


Equations  (C7)  and  (C8)  intersect  (in  the  7,  i  plane)  when 


—  =  1  -0-1(1  -2a)  (C9a) 

H 


Thus,  when  i  <  </2  the  buckling  load  is  given  by  Eq.  (C7)j  when  shear 
deformation  effects  are  no re  significant,  |  1/2,  the  buckling  load 

is  given  by  Eq.  (C8) .  The  huck'ing  load  >  is  plotted  as  a  function  of 
i  for  a  =  0 • 00549  in  Fig.  ly,  alhocgh  this  figure  is  applicable  as  long 
as  Q  <  <  1.  For  the  she', 5  comparab !e  to  that  described  by  Eqs.  (8.1)  and 
(8.2)  (i.e.,  §  =  0.1367)  the  buckling  load  is  given  by 

7"  «  0.837  (CIO) 

The  corresponding  critical  pressure  is  p  =  7078  k  10^  psi 
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APPENDIX  0.  PERTURBATION  SOLUTION 


The  result'-  described  in  section  8  show  that  the  stresses  are 
approximately  linear  functions  of  the  parameter  >  (see  Figs.  3  to  10).  This 
observation  suggests  that  a  perturbation  solution  should  give  accurate  pre¬ 
dictions  of  the  stresses  and  displacements  when  the  parameter  ^  is  small 
Isee  Eqs.  (4.5d),  (2.  I  )  and  (2.2  )}  For  the  initial  deflections  con¬ 
sidered  here  [see  Eq-  (4.1)],  the  governing  differential  equations  [Eqs. 

(2.  1  )  and  (2.2  )]  are 


w,xxxx  +  (VB)(0/L)hw  =  (4/B)(e/L)‘twp  -  4(0/L)2 


XX 


a)  ^  =Bw  -  4£(0/L)+(w  -  wp) 


(D  la) 

(D  lb) 


where  wp  may  be  obtained  from  Eqs.  (4.2). 

When  •}•  is  snisll,  a  solution  may  be  assumed  to  be  of  the  form 


w  =  L  w.  (x)  ' 


i  =0 


u)  =  Z  a),  (x)ili 
i  =0  ' 


(D2) 


Substitution  of  Eqs.  (D2)  into  (Dl)  leaves 


I"o,xxxx  +  (e/04wo  -  <V8)(8/L)\]  ♦  +  MKwoL 


+  (4//B)(0/L)2wo  ]  +  *2[  ]+•*•-  o 


(D3a) 


and 
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l®. 


C,K 


•  fiw 


*  44(0/1)  (w  -w  )]  +  *[®.  -Bw.  +  4| (6/ L)  w.; 

Op  I  f  X  I  y  XX  I 


+  i'2! 


=  c 


(D3b) 


The  first  approximate  solution  wq,  oj^  may  be  obta’ned  by  setting 
i/  =  0  in  Eqs .  (03)  •  Thus, 

w  +  (4/B)(0/L)4w  =  (4/B)  (0/L)4w  (D4b) 

o,xxxx  v,'o'p 

03  =Bw  -  45(8/L)i+(w  -w  )  (D4b) 

OjX  o,xx  '  o  p 

The  solution  for  wq  is  taken  in  the  form 

"o  ■  V  1  *  CI0?|W  *  c2ohM  *  c30»}<x)  +  CW9„W) 

where 

g  I  (x)  =  sinh(0x/4/B)cos(0x/t/B) 

92  (x)  =  s  inh  (0x/k/B)s  in  (Sx/L/B) 
g^(x)  *  cosh  (0x/L/B)  s  in  (0x/l/"B) 


(05) 


g^(x)  *  cosh (0x/l/B) cos (0x/l/B) 


I w  • 


The  arbitrary  constants  c.q  are  determined  by  using  symmetry 
[see  Eq.  (4-9)]  and  by  applying  Eqs.  (2.9b,c)  in  conjunction  with  Eqs. 
(3- 18),  (3.20),  (A. 34),  (A. 35),  (4.10),  (D4)  and  (05).  Thus, 

t S ( /ft 11 J i  t  -  [  (7i-ra,)/(ra-Mp) 

7  t?j(x)/B)  (VV’O 


and 


u;o  =  -  (32A*/L2)(7A-7A,)l-[0/2B)  +  t  (G^/G^)  H  92  (x)/B] 

3X 

+  [  (I/2B)  (G40/G20)  -  Ug4(x))  (d6) 


who  re 


G10  (2/q)  (g^0‘*'9j0)/t  B;  ^^920  +  9)0%0)  +  B2 (g30940  "  910920^ 

G20  =  ^8I9I0  +  B2S50^/'-Bl  ^g30920  +  9I0940^  +  B2  ^930940  "  910920^ 

G4Q  =  (B>930  "  B29!C)/lB!(930920  +  9iG940)  +  B2(g3o94o  '  910920}j 
A*  =  ^0G2O/{^*  +  GJ0(l-ec)  tp"  +  (l+6-26Vr)ac]/3cl 

g.  =  g. (L/2) 

3io 


(continued  on  next  page) 
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6,  =  I  +  2£(!  -  25  +  a) 

b2  =  1  -  2|(l  +  2i  -  a) 

The  corresponding  stresses  may  be  obtained  by  direct  substitution 
of  Eqs .  (D6)  into  Eqs.  (53).  This  solution  could,  of  course,  be  obtained 
fro.n  the  solution  presented  in  the  text  by  simply  setting  i|i  =  0  in  Eqs. 
(4.14)  and  (4. 13) • 

The  first  perturbation  yields 

+  (VB)(e/L)4t,  -  -  (4//»H0/o\„  (»?•) 

and 

“i,»  '  “i,„  -  »>7b) 

with  boundary  conditions 

w,  (±  L/2)  =  0  ,  »  (±  L/2)  =  0  (D8) 

The  solution  for  Wj  can  be  shown  to  be 

w,  =  -[8(7A->'A1)/e2)(c21g2(x)  +  c^g^x)  +  (0x/L)  [  c ,,  g  (  (x)  +  Cj(9j(x)]J 


(09a) 


where 

CI1  =  "  *  ”  ^Gi4(/G2Qy  J  (D9B) 

Cjj  =  -  0/2)1  I  *  CG4o/C2o) j  (09c) 

C21  =  *  G40  +  B  WlBl  (%o22C+9|0^0)+B2(  Vu0-9|0920}  j  (D9d) 

C4l  =  +  *  G20~B  920/[B1(9^20+^I0940)  +  B2  bo%Q~*  I  O^O3  J  (D9c) 

A  =  (-6/2)  (c  j  j9  j0  +  c^g^)  (d9<) 

B  =  -  [  l-26(2S-a)jlc),g10+c5)g50+(e/2)lc,  ,  (o^-g^)^,  (920+94o)  0 

-  2^"2c|  i930+2c5,g]0+^/?)i-c) ,  (920+940)+c^,  (94o-920)  J]  (D9g) 

The  solution  may  now  be  written  to  include  the  first  perturbation, 

thus 


=  w  +  itw. 
o  I 


w 


and 


Ul 


(DIO) 
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When  y  and  £  are  small  (as  they  are  in  the  region  of  p  for  which 
the  shells  considered  in  this  work  behave  elastically) 


B  =  1  -  4/  4  ~  * 


[see  Eqs.  (2.4)) 


Then  the  linearized  form  of  Gjq,  G^q  and  become 
GI0=FI1'  '«  1 

C20  =  f2  -  »  (S  *  f2  F3  > 

G40  -  F4  +  2^F2  -  3^  FJ  F4  J  (D>  •) 

where 

F|  =  +  9]o^^3o920+91094C^  +  ^950940"910920^ 

F2  =  ^910+950^/^930920+s1C940^  +  ^930940"910920^ J 

0^F3  =  [940930‘910920  *  %0920+910940)]/[(930920+910940)  +  (930940-910920)] 
\]  3‘4i"  F4  =  ^930‘910^^930920+910940^  +  ^93094o"910920^  ^D'2^ 

The  functions  F.>  correspond  to  the  functions  F,  listed  in  Ref- 
(12)  which  were  displayed  graphically  as  functions  of  0  for  the  isotropic  case. 
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The  stresses  can  be  obtained  by  substituting  Eqs .  (DIO)  into 
Eqs .  (53)  arid  using  Eqs.  (5-6).  In  general  Eqs-  (Dll)  and  (D 1 2 )  may  be 
used  in  place  of  the  definitions  for  G  q  given  in  Eqs-  (D6) .  When  this 
is  done,  the  charts  of  F.  v<;  9  offered  in  Ref.  13  may  be  utilized  to 
obtain  G  This  procedure  results  in  a  considerable  reduction  of  compu¬ 
tation  time  and  provides  accurate  stresses  and  displacements  for  the  most 
general  case  provided  that  i)i  is  reasonably  small. 

To  obtain  the  F.  vs  9  relations  displayed  as  functions  of  y  as 
well  as  9  in  Ref.  6,  it  is  necessary  to  set  (  -  6  =  a  -  0  in  Eqs.  (4.12)  and 
to  replace  Gj,  G^  G^  and  in  Eqs  .(4-12),  wi  th  F  j ,  F2,  \J  0-91/3  F^  and 
J^0.91/5  F^,  respectively.  These  functions  are  displayed  graphically  in 
Ref  5. 
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APPEND!!'.  E 

STRESSES  IN  NONHQMOGENEOUS  SHELL 

It  was  observed  in  Ref.  that  because  of  the  nonhomogeneity 

of  material  properties  through  the  thickness  of  the  composite,  Mill's  yield 
criteria  for  a  homogeneous  orthotropic  material  could  not  be  expected  to 
provide  accurate  estimates  of  the  yield  pressure.  This  short  coming  moti¬ 
vated  the  work  that  is  discussed  in  this  appendix. 

The  elastic  constants  given  in  Eq .  (8.1)  were  obtained  experi¬ 
mentally  by  measuring  the  strains  in  a  pressurized  cylinder  (of  JC:2L  winding 
ratio)  and  utilizing  stress  strain  relationships  in  a  simple  homogeneous 
shell  theory  (Ref.  14).  Hence,  it  is  reasonable  to  expect  that  the  homo¬ 
geneous  shell  theory  utilized  in  this  work  will  accurately  predict  the 
shell  deformations  [Eqs-  (4.13)  and  (4.14)].  However,  the  nominal  stresses 
obtained  with  these  constants  for  an  equivalent  homogeneous  material  [Eqs. 
(5-3)]  cannot  always  be  expected  to  correspond  to  the  actual  stresses  acting 
in  the  nonhomogeneous  composite  shell,  which  consists  of  materials  of  widely 
varying  elastic  properties. 

If  the  deformations  .redicted  by  shell  theory  are  accurate,  then 
the  individual  constituents  of  the  shell  material  must  deform  according  to 
Eqs.  (4.13)  and  (4.14).  Once  these  deformations  are  known,  it  is  necessary 
to  use  Eqs.  (5-0  to  calculate  the  stresses  in  the  constituent  materials- 
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It  is  possible,  however,  to  obtain  approximate  values  for  the 
stresses  in  the  nonnomogeneous  shell  by  utilizing  the  uniaxial  stress 
strain  relationships  suggested  in  Ref-  15  and  Ref-  |6.  For  example,  the 
effective  composite  modulus  Sf'V>  *n  a  direction  normal  to  the  fiber 
direction  could  be  obtained  by  assuming  the  fiber  resin  system  to  behave 
as  a  spring  in  series.  This  results  in  (see  Ref-  15) 


1+V  /V 


^  £  / — - 3. — £ — 

c2  e  lt£  V  /E  V 
eg  g  e 


)  for  >  68  percent 


(E.l) 


in  which  E^  and  E^  are  the  moduli  of  the  resin  and  fiberglass,  respectively, 
Ve  and  Vg  are  the  percentages,  by  volume,  of  the  resin  and  fiberglass.  In 
Ref.  15  it  is  shown  that  a  mere  complicated  relation  is  necessary  for  com¬ 
posites  with  Vg  «r  68  percent. 

It  is  noted  that  the  effective  modulus  given  by  Eq .  (E.l)  is 
based  upon  the  assumption  that  the  spacing  between  the  fibers  is  uniform  in 
that  the  volume  percentage  ratio  v^/V^  is  assumed  to  be  equal  to  the 
(presumably)  uniform  ratio  of  fiber  cross  section  to  resin  cross  section. 
Shaffer  has  pointed  out  winding  irregularities  could  cause  serious 
fluctuations  in  this  ratio. 

The  effective  composite  modulus  Ej  in  a  direction  parallel  to  the 
fiber  direction  is  given  by 


1+E  V  /E  V 

E  =  E  ( - e-  £  -3-3.) 

t  j  \i  /u  j-i  > 


V  /V  +1 
e  g 


(E-2) 


Shaffer,  B.W.j  Filament  Reinforced  Plastics  Micro-Mechanics-Structural 
Mechanics:  Seminar  given  at  Polytechnic  Institute  of  Brooklyn,  Feb.  9,  1 9&S . 
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!n  the  fiiament  winding  direction  the  stresses  ineach  of  the  con¬ 
stituents  may  be  obtained  by  multiplying  the  stresses  in  the  equivalent 
homogeneous  shell  by  appropriate  ratios  of  the  elastic  constants  (these 
ratios  are  often  referred  to  as  stress  concentration  factors).  There  is 
some  approximation  involved  due  to  neglecting  Poisson  ratio  changes,  but 
the  error  induced  in  this  manner  is  expected  to  be  small.  If  more  accuracy 
is  desired,  the  deformations  given  by  shell  theory  could  be  substituted 
directly  into  Eqs.  (5-0-  In  the  directions  normal  to  the  filament  winding 
direction,  the  stresses  are  transmitted  from  one  constituent  to  the  other 
and  must  be  matched  across  the  glass  resin  interface  in  order  to  satisfy 
equilibrium  requirements.  In  these  directions  it  may  be  assumed  that  the 
stresses  in  the  fiberglass  and  resin  are  equal. 

Since  the  stresses  in  an  equivalent  homogeneous  shell  have  been 
calculated  (see  Table  3),  the  stresses  of  the  constituents  may  be  approxi¬ 
mated  as  fol lows: 


1.  In  a  fiberglass  filament 


ah  E  ( 1  -V  )  ah 

(.24  c  -S - -  (-*_) 

Jg  c  2.  VN  ’ 

*  E  (l-y  )  x 


oh  1+V  /V  E  oh 

}g  =  <“£  V  /t~]t 
x3  eggecpx 


in  a  longitudinal  layer 


(continued  on  next  page) 


84. 


oh  i+V  /V  E  o  h 

Vi  ;g  N+E  V  /E  V  ;E  S  ' 
x  3  e  g  g  e  x  x 


in  a  circumferentially  wound  layer 


<U>  E„(l-v  )  ah 


<s** 


(ji^) 


V-.  >  * 


(E.3) 


2.  tn  the  resin 


oh  E  ( I  -v  )  a  h 

v  a  '  ’  v 


<fT-) 


e  "  2,  ^ 

£x(l-vt  )  x 

A  t 


ah  1+V^/V  E  oh 

*N  Je  =  (1+E 

x  e  g  g  e  <p  x 


a  h  1+V  /V  E  o  h 

=  / - 9__£ - v_e  ,_x_) 

'N  ;e  '•l+E  V  /E  V  'E  'N  } 
x  e  g  g  e  x  x 


in  a  longitudinal  layer 


in  a  circumferentially  wound  layer 


oh  E  .  2  oh 

<N^e  ■  E* 

X  <p  1  -v  X 


(E.4) 


V  and  ve  are  the  Poisson  ratio's  of  the  giass  anc  resin-  The  g  and  e  sub¬ 
scripts  denote  the  local  stresses  in  the  fiberglass  and  the  resin,  respectively. 

Once  the  stresses  have  been  obtained  from  the  above  relations  a 
maximum  stress  yield  criterion  can  be  used  to  approximate  the  magnitude  of 
the  pressure  at  which  yielding  begins-  Thus, 
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2(h/a) 


(l+6)2(“) 

X  5 


in  the  fiberglass 


Py 

=  2(h/a) 

Ye  (W)2^)a 


in  the  resin  (E.5) 


where  Yg  and  Ye  are  the  yield  stresses  cf  the  fiberglass  and  the  resin., 
respect i vely . 

Limited  calculations  were  performed  on  o  shell  of  dimensions  equal 
to  those  of  the  equivalent  homogeneous  shell  [Eqs.  (8.1)].  The  constituent 
elastic  constants  and  yield  stresses  were  taken  from  Ref.  (18),  while  the 
volume  ratio  was  deduced  from  Ref.  1. 

Eg  =  10-5  x  106psi  ,  Ee  =  0.5  x  lO^psi 

vg  =  0.2  ,  ve  =  0.35 

Y  >  200,000  psi  ,  Yc  =  22,000  psi 

V  /V  =2. 188 
9  e 

* 

Substitution  of  Eqs.  (E .6)  into  (E-5)  to  (E-3)  gives  approximate 
values  of  the  initial  yieid  pressure.  It  is  found  that  the  yielding  would 
initiate  in  the  resin  near  the  shell  inner  surface,  next  to  the  frame,  at 
pressures  of 
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Py  =  11785  due  to  circumferential  stress  in  a  longitudinal 

e  layer 

Py  =  11157  due  to  longitudinal  stress  in  a  circumferential 

e  layer 

These  yield  pressures  are  within  ten  percent  of  the  experimentally 
determined  collapse  load  of  the  shell  examined  in  the  text-  Hence,  the 
shell  appears  to  be  weakest  at  the  frame,  near  the  inner  shell  surface. 

The  stresses  in  the  direction  perpendicular  to  the  fiber  directions  exceed 
the  relatively  low  yield  stress  of  the  resin  at  extremely  low  pressures- 
When  this  happens  the  shell  might  be  considered  to  be  delaminating  in  that 
the  resin  is  no  longer  binding  the  layers  together. 

Winding  irregularities,  as  noted  earlier,  may,  through  effective 
modulus  changes,  cause  yielding  to  occur  very  early  in  some  portions  of 
the  shell.  Evidently,  this  early  yielding  results  in  little  more  than  slight 
changes  in  the  composite  moduli  (e.g./  see  Ref.  14  and  Ref.  19). 


TABLE  ! 

STRESSES,  DEFLECTIONS  AND  LOADS  (m  =  5  =  P  -  0) 
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TABLE  2  Vi  ELD  PRESSURE 


- r 

YIELD 

CRITERION 

P  -  PLASTIC 
POTENTIAL 

YIELD  STRESSES 

X  io"5  PSI 

Y!  EL 

PRESSURE  x  10" 

4  PSI 

1 

Mi  DBA Y 

OUTER 

SURFACE 

Mi  DBA Y 
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SURFACE 

M-  MAXIMUM 
STRESS 

Y 

X 

Y 
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B 

h. 

i  n . 

_ 1 
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0.194 
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0-388 

ORTHOTROPIC: 

E 

X 

=  4.74 

x  10 

A 

PSI,  E 

=  6.  i4 
s 

x  106  PSI,  Ef  = 
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o 

it 

CL. 

1.0 

1-7 
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0.3992 

WBPPVI 

0.4325 
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o 

II 

N 

t> 
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a. 
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OO 
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0-37*0 
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0-3979 

■BBB 

■'ML  I'll 
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0 . 2325 

P(az  i  0) 

.T’ 

.C 

0.22 

0 . 3740 

0. 1753 

0.4332 

0.4997 

0  -2503 

M 

1.0 

1.7 

1.858 

0.6641 

2.637 

1 

i 

ISOTROPIC: 

E  - 

E  -  E 

X 

=  5-355  x  106  PSI.. 

Bf  --  1.1 

6 

x  iO  PSI 

! 

P  (°2  =  0) 

1.0 

1.0 

0.22 

0.3753 

0. 1700 

0.4015 

0.1921 

0-4373 

0-2397  j 

P  (oz  =  o) 

1-7 

1-7 

0.22 

0-3754 

0. 1705 

0 . 40 1 5 

c. 1922 

0.4383 

C.24CU  j 

P(az  =  0) 

i} 

0.22 

0-3734 

0. <707 

0.4021 

0.1922 

0.4389 

0.240 ;  I 

P(a2  *  0) 

O' 

1 

0.22 

0.3754 

c. 1707 

0.5060 

0 . 2  i  15 

0 . 5054 

0.2533  ! 

P (az  =  0) 

1.0 

1.0 

1.0 

1.655 

0.7157 

1.771 

0.8718 

1.855 

1 . 05,.  | 

P  (az  =  0) 

1.7 

1-7 

1-7 

2.882 

1.217 

3.01 1 

1 .482 

3-153 

l  785  i 

f 

M 

1.0 

]  .0 

1.808 

0 . 6503 

1 .610 

0-9004 

l 

M 

1-7 

1-7 

3.074 

1.071 

2.756 

1-551 

1 

1 

Maximum  Stress  -  Axial 


Maximum  Stress  -  0 i rcunferen L i a i 


Reproduced  From 
Best  Available  Copy 
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table  2  YIELD  PRESSURE 


YIELD 

CRITERION 

P  -  PLASTIC 
POTENTIAL 

M  -  MAXIMUM 
STRESS 

YIELD  STRESSES 

X  lO-5  PSI 

YIELD  PRESSURE  x  10‘ 

4  PSI 
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OUTER 

SURFACE 

MIDBAY 

MEDIAN 
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H 

n 

Y 

.z 

_ 

h, 
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l 

B 

=  Ef  =  5-395  x  !06  PSI 

P(o2  =  0)  ' 

1.0 

1 .0 
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0-3946 
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0 . 2304 

P(o2  =  0) 

1-7 

1-7 

0-22 
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0.3953 
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0.4242 
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1.0 
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1-7 
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2.98I 
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Z 

V; 

ft 

-O' 

M 
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0.6825 
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M 
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■ 
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TABLE  3- 


TRANSVERSE  SHEAR  DEFORMATION  AND 


BEAM  COLUMN  EFFECTS 


(A  =  0,  62 
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0 
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i- 
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SCALE  SECTION  OF  STRUCTURE  (WITH 
RECTANGULAR  CROSS  SECTION  RINGS) 


bo 


Po 


J  I  •  I  J 


TT 

lO 

d  bwr 

liliiiini 


JLtU.ll  I  It 

tTTTI 


T 


TTT 

Po 


T 

h0 


-Pi' 

HI  iit 


bi 


«)  CROSS  SECTION 
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c)  FREE  BODY  DIAGRAM  OF  SECTION 
OF  SMELL  IN  CONTACT  WITH  RING 


FIG.  2  REINFORCING  RING 
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FIG.  3  AXIAL  STRESSES  AT  FRAME 


xIO  (PSD 


r 


FIG.  4  CIRCUMFERENTIAL  STRESSES  AT  MIDBAY 
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FIG.  6  CIRCUMFERENTIAL  STRESSES  AT  FRAME 


FIG.  0  AXIAL  BENDING  STRESSES  AT  FRAME 
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FIG,  9  AXIAL  BENDING  STRESSES  AT  MIDBAY 
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FIG.  to  CIRCUMFERENTIAL  BENDING  STRESSES  AT  MIDBAY 


FIG.  M  EFFECT  OF  INITIAL  DEFLECTION  ON  STRESSES 
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FIG.  12  EFFECT  OF  INITIAL  DEFLECTION  ON  YIELD  PRESSURE 


FIG.  13  EFFECT  OF 


CT  OF  INITIAL  DEFLECTION  ON  YIELD 


EFFECTS  OF  TRANSVERSE  SHEAR  DEFORMATIONS  ON  STRESSES 
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FIG.  17  EFFECTS  OF  INITIAL  DEFLECTION  ON  STRESSES 
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TABLE  3.  TRANSVERSE  SHEA A  l TFORMAT I  ON  AND  BEAM  COLUMN  EFFECTS 


V  =  0.2 

V  =  0 

MIDBAY  FRAME 

v  - 

MIDBAY 
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E  w 

X 

N 
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